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Department of Theoretical Physics Bolestaw Bierut University, Wroclaw 
(Received September 15, 1958) 


In this work an equation for multi — component lattice gas is obtained. The thermo- 
dynamical potential, the mean — value of energy, and the mean value of numbers of 
particles for high and low temperatures are calculated. The method used in this work is 


the variational method of Bogolubov with one variational parameter. 


1. Introduction 


In this work we consider the theory of multi-component lattice gas with the 
help of Bogolubov's variational method of computation of statistical sums. This 
method quoted in Kwasnikov's papers (1956, 1957a, 1957b) in the case of Ising's 
model of a ferromagnetic body (Kwasnikov, 1956) gives good results both for high 
and low temperatures. The remarkable feature of this method is the fact, that one 
have not to restrict him self to interactions vanishing at distances larger than those the 
nearest neighbours. 

The lattice gas model was introduced to creating an approximate theory of con- 
densed states (see Newell and Montroll, 1953; Yang and Lee, 1952). It is known that 
the statistical sum for a one — component lattice gas may be obtained by elementary 
calculations on‘ the basis of the statistical sum of Ising’s model of a ferromagnet.) 
(Newell and Montroll, 1953; Yang and Lee 1952). Holo un this method refers 
only to one — component lattice gas. 

The theory of the lattice gas was worked out by Levine, Mayer and Aroeste 
(1957) (one component gas), and Levine (1957), (multi — component gas). They 

obtained the first terms of an expansion in powers of 1/kT. Their methods proceeds 

in quite a different direction from that of Ising. This was resulted (i) from the fact 
that in Ising treatment there are no exact solutions for the three dimensional case 
and (ii) from the difficulties in introducing a simple approximate method. 

In this paper a description of the Ising's type is applied to the multi — component 

lattice gas with the use of Bogolubov's method of approximation. 
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2. Statement of the problem 


We consider an f component lattice gas. The symbols (7) ynd (j,,) denote the 
positions of the i-th particle of the J-th component, and j-th particle of the m-th 
component of lattice gas. The coordinate (ij) can take any one of Q discrete values. 
In this work the indices run as follows: À, i,j = 1, ... N, s, t, u = 1, ... Q k, l, m, = 
=P 

We define: 


Fx... Np = (3)... >). (3. Pl exp (— BE) 


(11) (Ny y) (Np) 


where summation over (ù) denotes summation over all possible positions, and 


x M 3 Ni Vee (ie hr) — : 2E > x Ad wa” (ir, hm) 


d h>i h,i 


where summation over d denotes summation over all possible distances in the lattice. 


Functions vår” (i, h„) are defined by: 


hal EL —oo when (i,) and (h,) are on the same lattice site 
Ag" Va Gp ha) = He when (7, and (h,) are distanced on d. 
O otherwise 

We define the value: 

E=-X3 D, AGG— > NCE (2.1) 

k (s>t);d k+m (s> t; d 

where: summation over (s > t) ; d denotes the summation over all possible pairs of 
lattice points separated by a distance d, and summation over all possible d. CË can 
take on either of two values: 0, 1. The case C* — 1 denotes, that at lattice point s we 
have an atom of the k-th component, the case C* — 0 denotes, that at lattice point s 
there are no atoms of the k-th component. We see, that E — E' for such a choice the 


variables C^ which corresponding to the actual configuration of the atoms if E oo. 
There are hovever two restrictions: (7) 


2) DE 


s km 
and (ii) 


2 G= 


The restriction (i) corresponds to the fact that each lattice point s is occupied by zero 
or one of the atoms; restriction (ii) corresponds to the fact, that in this lattice there 
are N, atoms of the k-th component. We have: 


7 Q-3-Qpecm- 


DD mm ap 


Sent ETIN AD) (yey yas a (>) - >) exp (BE?) 
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Å : x é 2 š Å 
where: summation over C? denotes summation over all possible combinations of 0 and 
l with restrictions (i) ynd (i). The grand — partition function is given by: 


NES 
ZN 
på me Ny... Nf 
N,= 0 Nf = l 
ZNM=2 
k 
EDIE 1. D 
N=0 Nf=0 m cl cf 
1 Q Q 
«ep (83 > HOW NEA GEOP) 
k (s>t);d k= m (s> t): d 
ADD) (02 202 PD Toe 
HE us ie 


pg GODE Ro OC) 


k (s7D;d R+ m (s> 1);d 
where: på = BAR; où" — BAR" and Z, = e^ 
It has been proved by Bogolubov, that for sums of the type: 
Spexp(— E) we have: 
Z = Sp exp (— E) > Zi = Z, €? 


where: 
1 7 F „ 
Zo—Spep(—Rh  S--.BS(E'ep(—R) Ë= K+E 
0 


In the method of computations used here we define E, and E" with help of the 
variational parameters. By means of differentiation with respect to the variational 
parameters we maximalize Z;,. The approximation consists in considering the ex- 
tremum Zi, as exactly equal to Z. 


3. Details of the calculations 


We choose | 
| ES c (3.1) 
|» s 


where: ©, = (uj — Ë) for one — parameter “simple“ theory, and ©, = (uj — 63) 
for f parameter theory. It is clear that the f parameter theory gives a better approxima- 
tion than the “simple“ theory. By virtue of (3.1) and (2.1) 
ree puro. (eme TD va; CP 
(s>t);d k#m (s> t); d 


£ = E, + E; + Es 3.2) 
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where: », = for “one parameter theory“ and v, = ¢, for “f parameters theory" 
Z= (>) sy ODIUM 
c! cf cl c s k 
(without restriction 71) 
=I [I - 2) exp (22 0] = D + >; e@%]° 
s ci cf k k 
(without A it) 
Z, S= Sp UU exp (23 E) — 4 + 4, d Az 
where: Ag = Sp {E’g exp (— E,)} 
We compute the values 4,, Ay A3 
Ay = (DJ ++ Dy) ++ (22 22 HT exp (>; 6, C) (— 22 22 G.) 
Gb cf pr Che k "p 
(without 77) 
--9»»(2.*2)Ge»(2;0, Cr} (1 + Dy Em)?! 
l cl cf k m 
(without ii) 


= — Q > vy, e& (1 + > Pod å 
4--—2 - 2 - 3X TAC Me Tc 
1 f ci ch d u k 


1 (sts 
Ci Ci Q 


(3.3) 


(without ii) 


=—92 2 va ra (X e D) (D 2) G G I exp (22 6, C2) 


l d 1 
š Ci 1 Ca Co 


er — Q D>; g! e(l + >) e68)-? 
7 k 


where: Qn; is a number of pairs of lattice points separated by an distance d, and 
= n 
yis 2 Pa "a 


In a similar manner we obtain 


As = — Q Y på Pkt Om (Len)? 
km l 


1 1 
Pint T DQ In Zune = ro In Z, — va) S = In U + eva (uv, + gl e8i) edı + 
+ U-2 >) pm e9k + Om e (3.4) 
km 
where U = (1 + >) e&) 
k - 
We now consider the “one parameter theory“. In this case: 


Vint = In (1 + ve) + t» ( E»)? He (ée 2 — (3.5) 
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where: v 2; ec 2i a X gna, 
k 


k=+ m 
We maximalize the yt. The condition (vi, PC) = 0 is equivalent to 


Zen 
Eee tn 3.6 
d yety CBS) 
When a ‘quadratic form & (in variables z,) is positive — definite, then £ is negative. 
With the help of (3.6) we cau express the maximum value of Yigg, by means of the 


parameter ¢ in a simple form. We have: 


e pte inm oA inf 4 (3.7) 


Wint = In [et (et + ») 


s 
222 Ae Ag DE 


We find here the average — number of particles of the /-th component. 
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In virtue of (3.8) 


N= YR = 2 69 
l 


We calculate now the average energy of lattice gas. 


a Gp). Te, n 
TOES T E 


In the case of a one — component lattice gas the “simple“ theory gives the same 
equations as the “f — parameter theory“. In this case: € = gz?: v = z 
Gå oe: Qc 


u ler ne 


z- opis (5) (3.11) 


4. Solutions of equation (3.6) 


We have found the approximate solutions of equation (3.6) in two limiting cases: 
i) (eh?) = (a/v) > 1 — the case of low temperatures 
ii) (2e/v?) = (a/v) > 1 — the case of high temperatures 
` From the form of equation (3.6) it follows, that this equation possesses only one so- 
lution for v > 0, « > 0. 
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(This is seen immediately by considering a plot of the function — œ (e^ + »)71) 


Let us consider the case i. We rewrite equation (3.6) in the form: 


We have: 


a a 
myr): Gic den 


We define the sequence C, by the formulae: 


& & 
ea TR et, 
We can easily prove that: 
)6«& i) huic 
Similarly, we have 
& = = & 
year and: Cn41 ce ES 


dues a eem MET eS 
We may evaluate the difference: | Gn 6,41] — 


ein — ein 


(en + v) (etn + v) 


M [64 & 
ein-F» dey 


— & 


1 = 1e es 
= (Eni Cn+1] len — ešn| < = (Cati Cnil ein = Cn+1 = 1 —< an 
By the above formula: 


aia rni leone tt 


k=1 
motn å Få OE 
= léo — Col I] E Ce lr Eu | (4.1) 
k=1 


It follows from i), ii) and iii) that the sequences €, and Š, converge. The sequence €, 
converges particularly rapidly if (x/v) > 1. In this case we can write B 


———  —— (a! 


(4.2) 
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In other cases we can use another approximation: 


== 
- Í = å 
e =s Ure eS 9 exp ap (—2) +5 +” si 
=y —1 
+ "Ces + y (4.3) 
Now, we consider the case II, in which (a/v) < 1 
Hetty =—(2)y=— on 0<1 (4.4) 


La E = > C, O”. Thus, we have: 


n=1 
EI 1 
er 6rn-6e (at pet) ot (nosse te) ot 
1 LM 165 
+ 04 (Ct at at art Ast) +00 


Puthing these expressions into (4.4) and equating the coefficients of equal 
powers of ©* of both sides we get: 


y : —y2 
NS Se en: 
253 iL y? 


—6 12631 


Č; = 


1 
Ra en 4 
— 3y4 zr t3 i v3 Lee 


"CO GE pH 3641) 


For one- component lattice gas we have in case i): 


for A>Oifk4T=0 N=Q 
I would like to thank Dr Lopuszanski for helpful discussions. 
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DIE ELEKTRONENMIKROSKOPISCHE UNTERSUCHUNG VON 
SPALTFLACHEN DER InSb-EINKRISTALLE 


A. FeLTYNOWSKI UND L. GORSKI 
Institut får Physik der Polnischen Akademie der Wissenschaften, Warszawa 


(Eingegangen am 15 November 1958) 


InSb-Einkristalle wurden långs der (111) und (100) Flåchen gespalten. Mittels eines 
geeigneten Abdruckverfahrens (Zaponlack-Aluminium) wurde die Struktur der Spalt- 
flåchen untersucht. Die Spaltstufen (cleavage steps) und die Gleitlinien sind die Hauptele- 
mente der Struktur der Spaltfläche. Die Spaltstufen entstehen infolge von Schrauben- 
versetzungen, deren Verschiebungsvektor senkrecht zur Spaltfläche steht, bzw. in allge- 
meinem Fall die Spaltfláche schneidet. An der Korngrenze bilden sich meistens kleine 
Spaltstufen, die sich später zu grösseren Stufen vereinigen, und sogenannte „Flussfigu- 
ren“ (river patterns). Die Spaltstufen verlaufen in Richtung der Spaltverbreiterung. 
Auf der Spaltfláche (100) beobachtet man zwei Gleitliniensysteme (mit dem gegenseitigen 
Winkel 45°), die von zwei verschiedenen Gleitflächen (111) und (110) herrühren. 


Die vorliegende Arbeit ist mit den Untersuchungen über verschiedene Eigen- 
schaften von Einkristallen der wichtigen Halbleiterstoffen verbunden, was ein Thema 
der komplexen Arbeiten unseres Institutes darstellt. 

Das Ziel dieser Arbeit war festzustellen, wie weit die Methode der elektronen- 
mikroskopischen Untersuchung von Spaltflächen eines Einkristalls zur Untersuchung 
der Strukturdefekte mit deren Einfluss auf Eigenschaften der Halbleiter geeignet ist. 

Natürlich die Sichtbarmachung der Defekte durch die Âtzung der Oberflüche 
der Einkristallen und nachfolgende lichtmikroskopische Beobachtung ist seit einigen 
Jahren bekannt und in sich eine einfache Methode darstellt. Vogel (1953) und Ellis 
(1955) deuten die verschiedenen beobachteten Atzfiguren bei Germanium mit den 
Schrauben und Kantenversetzungen, Vogel und Lovell (1956) hatten mit dieser 
Methode Silizium untersucht; Indium-Antimonid ist letztens von Bardsley und 
Bell (1956) untersucht worden. | 

Die Atzfiguren scheinen keine submikroskopische Struktur zu haben und deshalb 
liegen keine elektronenmikroskopische Untersuchungen vor. Zwar ist hier eine Arbeit 
von Allen und Smith (1956) zu nennen, diese Autoren haben aber Rastermikroskop 
benutzt, Abdruckverfahren ist nicht verwendet worden. 
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Wir haben die Ätzfiguren sowie beim Germanium, wie beim InSb mittels Ab- 
druckverfahren elektronenmikroskopisch untersucht, die Resultate sind aber nicht 
interessant, da die Ätzfiguren keine Substruktur aufweisen. 

Deshalb scheinte uns sehr interessant die Spaltfläche der Halbleitereinkristalle 
zu untersuchen. 

Bei Germanium und Silizium kann man eigentlich nicht von der Spaltung 
sprechen, die Bindungsenergie ist zu gross und entsteht bei den Versuchen ein 
Durchbruch statt der Spaltung. InSb-Einkristalle lassen sich in den Flächen (100) 
und (111) leicht zu spalten, die Qualität des Einkristalls spielt hier auch eine 
wesentliche Rolle. 

Die Spaltfláchen von Einkristallen sind, so wie uns bekannt ist elektronenmikro- 
skopisch bis jetzt noch nicht untersucht worden. Die lichtmikroskopischen Beobachtun- 
gen an Spaltflächen betrafen lediglich der Einkristalle von Metallen, hauptsächlich von 
Zink, wo die Arbeit von Gilman (1956) zu nennen ist. Nach Gilman infolge der im 
Einkristall bestehenden Defekte bilden sich bei der Spaltung auf der Spaltfläche die 
Stufen (cleavage steps), die wir als Spaltstufen nennen wollen. Bilby und Smith (1956) 
bei der Untersuchung der Gleitebenen in den gespalteten Zinkeinkristallen nach deren 
Verformung haben so genannte Flussfiguren (river patterns) beobachtet, d.h. grosse 
Anzahl von Spaltstufen an der Korngrenze. 

Die von uns untersuchten Einkristalle von Indium Antimonid sind nach der Me- 
thode der Zonenschmelzung erzeugt worden. Nach der Feststellung der Kristallorien- 
tierung sind sie in verschiedenen Kristallflächen gespaltet worden. Zu diesem Zweck 
ist der Kristall in einer Goniometerfassung zu befestigen und mit einem Vidiamesser 
zu spalten. Zu den nachfolgenden licht — und elektronenmikroskopischen Untersu- 
chungen eigneten sich die Kristalle, die in den Kristallfláchen (111) und (100) ge- 
spaltet worden sind, wo zu bemerken ist, dass die Spaltung gemäss der Fläche (100) 
leichter durchzuführen ist. 

Die Spaltfläche ist erst unter dem metallographischen Auflichtmikroskop durch- 
zusuchen. Die Abbildung 1 zeigt beispielweise die reiche Struktur eines Teils der 
Spaltfläche bei der Vergrösserung 540 x. Aus der so untersuchten Oberfläche, die 
interessante Struktur aufweisen scheint, ist ein Abdruck zu machen, um diesen Teil 
der Spaltfläche im Elektronenmikroskop durchsuchen zu können. Wir haben leider 
kein Rückstrahlungselektronenmikroskop und mussten das Abdruckverfahren be- 
nutzen. Es war ein zweistufiges (positives) Saponlack-Aluminium-Abdruckverfahren 
verwendet worden (Saponlack von der Spaltfläche abgerissen, nachher mit Chromium 
unter 15° schattiert, nachfolgend mit Aluminium unter 90° bedampft; der Saponlack- 
film im Amylacetat gelöst, es blieb als Präparat nur mit Chromium schattierte Alu- 
miniumfolie). In manchen Fällen ist ein Mikromanipulator benutzt worden. Die in- 
teressanten Teile des Abdruckes sind mittels der Glassonde oder einer Nadel noch 
bevor der Trennung des Abdruckes vom Kristall zu zeichnen. 

Auf der Spaltfläche (100) des untersuchten Kristalls sind deutlich die Spaltstufen, 
. die in den Defektstellen der Kristallstruktur entstehen, zu sehen. Die Ursache dieser 
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Verankerung der Spaltung und nachfolgend der Entstehung einer Spaltstufe ist eine 
Schraubenversetzung, deren Verschiebungsvektor (Rurgers Vektor) sehnkrecht zur 
Spaltfläche steht, bzw. in allgemeinem Fall die Spaltflåche schneidet. Die Spaltstufen 
verlaufen in der Richtung der Spaltverbreiterung, wie auch Gilman in seiner Arbeit 
(1956) daraufweist. Ausser diesen Teilen der Spaltflåche, die mit den Spaltstufen 
verhältnismässig dicht bedeckt sind, sind grosse Teile (etwa ?/, der durchgesuchten 
Oberfläche) der Spaltfläche vorhanden, die ganz glatt sind und auf welchen nur sehr 
selten einzelne Spaltstufen zu finden sind. Diese Spaltstufen entstehen infolge der 
Verankerung der Spaltung durch nur einzelne Schauberversetzungen. 

Die ersten Bilder zeigen die Spaltstufen, die über ziemlich grosses Gebiet der 
Spaltfläche verlaufen (etwa 1004; Abb. 2 — 5). Die Spaltstufen verlaufen in der 
Richtung der Spaltverbreiterung. Am Anfang sind sie klein und dicht gelegen und 
nachher vereinen sich allmählich in immer gróssere, naturlich mit der abnehmender 
Dichte. Am Ende dieser Spaltstufen bildet sich auf der Spaltflåche ein grosser Struktur- 
defekt, der wieder ein Anfang anderer Serie der Spaltstufen ist. 

Auch die Korngrenzen sind oft ein Anfang oder eine Ende der Spaltstufenserie, 
wo die Spaltfláche auf anderen Seite der Grenze ganz glatt ist. Je grósser die Desorien- 
tierung der Körner ist (d.h. je grósser der Winkel zwischen den Körner ist), desto 
dichter die Spaltstufen an der Grenze. An der Korngrenze bilden sich meistens 
kleine Spaltstufen, die etwas weiter in gróssere sich vereinen, und bilden sogenann- 
te „Flussfiguren“ (river patterns). Deutlich zu sehen ist, dass die Dichte der Spaltstufen 
mit der Entfernung von der Korngrenze abnimmt (Abb. 6—8). 

Auf einigen Bilder (Abb. 9) ausser der gewóhnlichen Korngrenze ist auf der 
glatten Seite der Spaltflåche eine parallel zu dieser Grenze liegende Linie zu sehen mit 
den anfangenden Stufen, die nachher verschwinden. Das kónnte man mit der Tat- 
sache erklären, dass hier die Versetzungen des umgekehrten Zeichens hervortreten, 
die gegenseitig sich auslóschen. Solche „Auslöschungslinie“ kann man auch auf den 
glatten Teilen der Spaltfläche beobachten, wo sie mit der Korngrenze nicht in Verbin- 
dung steht (Abb. 9a und 9b). 

Auf einer der Bildserien hat man die Spaltstufen auch auf beiden Seiten der 
Korngrenze beobachtet. Am Anfang bilden sich die Spaltstufen selbständig (nicht auf 
der Grenze) und ziemlich dicht (Abb. 10) und finden ihr Ende: auf der Korngrenze 
(Abb. 11), wo die Dichte der Stufen wührend ihres Verlaufs stark abnimmt. Auf beiden 
Seiten der Korngrenze (Abb 11) sind die Stufen zu sehen, auf einer Seite helle, gróssere 
mir kleiner Dichte, auf anderer dagegen dunkle, kleinere mit grósserer Dichte. Dieser 
Unterschied, wegen der schrügen Bedampfung sichtbar, weist auf eine Desorientierung 
der Flächen auf beiden Seiten der Grenze. Es ist zu vermuten, dass die Höhe einer 
hellen Stufe der Hóhensumme der einzelnen dunklen Stufen, die zwischen den be- 
nachbarten hellen sich finden, gleich ist. 

In einigen Fållen konnte man die Änderung der Spaltstufenrichtung beobachten 
(Abb. 12—13). Die Ursache dieser Änderung ist ein Hindernis in Form von einer 
Versetzungsgrenze mit grosser Anzahl von Versetzungen, welchés die bisherige Rich- 
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tung der Spaltverbreiterung ändert. Gleichzeitig hat man auf dieser Grenze „Flussfi- 
guren* beobachtet (Abb. 14). 

Auf der Spaltfläche (111) hat man auch die kleinen und grossen Spaltstufen 
beobachtet und es war auch zu finden eine Serie der kleinen und dichten Spaltstufen, 
die auf einer Grenze entstehen. Muss man aber bemerken, dass die Spaltung entlang 
der (111) Fläche schwerer durchzuführen war, das Bild der Spaltfláche mehr zerrissen 
war, was mit Verbindung mit der niedrigeren Qualität des betreffenden Einkristalls 
weniger wertwolles Beobachtungsmaterial lieferte. 

Die lineare Dichte der Spaltstufen berechnet aus dem gesamten Bildmaterial 
beträgt von 104 bis 105 pro cm. Da die Anzahl der Spaltstufen der Zahl der Schrauben- 
versetzungen für durchgesuchtes Gebiet des Einkristalls entspricht, kann man die 
lineare Versetzungsdichte der untersuchten Einkristalle als von Gróssenordnung 10*— 
105/cm angeben. 

Auf einer Reihe von Bilder ist leicht zu bemerken, dass die einzelnen, grósseren 
Spaltstufen die lamellenfórmige Struktur aufweisen, die aus dichten und winzigen 
Linien besteht, welche oft unter dem rechten, oder nahe rechten Winkel zum Spaltstu- 
fenverlauf stehen. Es sind die Gleitlinien (Abb. 8, 9, 15, 16). Es sind auch die Gleitli- 


Ap». 1. Lichtmikroskopisches Bild von (100) Spaltfläche des InSb-Einkristalls. Vergr. 540x 


~ 


nien beobachtet worden, die unabhüngig von Spaltstufen verlaufen und bilden die 
Streifen von dichten Linien, die am meistens sehnkrecht zu Spaltstufen verlaufen. 
Diese Linien sind noch kleiner und dichter gelegen als die Spaltstufen. 

-In unserem Fall der Spaltung entlang der (100) Flüche entsteht die Gleitung 
hauptsächlich an den Flächen (111) und (110). Einen interessanten Fall stellt ein Bild 
(Abb. 16) vor, wo eine Spaltstufe von grösserem Ausmasse mit der winzigen Ultra- 
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Abb. 2 bis 5. Eine Spaltstufenserie, die über ein Gebiet etwa 1004 gross auf der (100) Spaltfläche des 
I  InSb-Einkristals verläuft. Vergr. 7.500 x 
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Abb. 6. Die Flussfiguren an der Korngrenze auf der (100) Spaltfläche des InSb-Einkristalls. Vergr. 
7.500 x 
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Abb. 7. Die Flussfiguren an der Korngrenze auf der (100) Spaltfláche des InSb-Einkristalls. Vergr. 
å 75500 xs ; 
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Abb. 9. Die Kornerenze und die „Auslöschungsgrenze“ Vergr. 7.500 x | 
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Abb. 9a u. 9b. Die einzelne Spaltstufe mit den Verankerungen und einer „Auslöschungsgrenze“ 
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Abb. 11. Die »Weissen* und "dunklen" Stufen auf beiden Seiten der Korngrenze. Vergr. 7.500x 
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Abb, 12 u. 13. Die Anderurlg des Spaltstufenverlaufs auf einer Versetzungsgrenze. Vergr. 7.500x 
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Abb. 14. Die Änderung des Spaltstufenverlaufs auf einer Versetzungsgrenze. An dieser Grenze sind 
. auch „Flussfiguren“ zu sehen. Vergr. 7.500 x 
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Abb. 15a. Die Gleitlinien und Spaltstufen auf der Spaltflüche. Vergr. 7.500 x 
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Abb. 16. Zwei Gleitliniensysteme auf der (100) Spaltflåche. Vergr. 7.500x | 
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Abb. 16a. Die Gleitlinien auf der Spaltflåche. 
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struktur der Gleitlinien zu sehen ist. Es sind auf diesem Bild zwei Gleitliniensysteme 
zu sehen mit dem gegenseitigen Winkel 45°. Das kónnte entsprechen dem Winkel 
zwischen den Schneidungsspur der Flüchen (100) (111) und der Schneidungsspur 
der Flüchen (100) (110). Da eine Gleitlinie die Schneidungsspur der Spaltfläche und 
der Gleitflåche darstellt, rühren die beobachteten zwei Gleitliniensysteme von zwei 
verschiedenen Gleitfláchen (111) und (110) her. 

Nach Mott und Nabarro (1948) liegt die Ursache der Gleitung in den Versetzungen, 
die im Einkristall schon bei seiner Erzeugung oder thermischen Bearbeitung entstan- 
den sind, wáhrend nach anderer Theorie von Kochendörfer (1938) die Ursache der Glei- 
tung in den Korngrenzen, oder in den Strukturdefekten anderer Art liegt, die bei 
kleinen Deformationsspannungen die Versetzungen erzeugen konnten. In unserem 
Fall liegt die Ursache der Gleitung wahrscheinlich in den Versetzungen, die infolge 
der Scherspannungen bei der Spaltung entstanden sind. 

Auf dem letzten Bild (Abb. 17) sehen wir einen Teil von grossen Graben (bzw. 
Höcker), der durch das betråchtliche Gebiet des Einkristalls låuft. Das ist wahrschein- 
lich ein inkludierter Kristallit, anders orientiert als der gespaltete Mutterkristall. 
Auf beiden Seiten dieses Grabens sind deutlich Gleitlinien zu sehen. 

Zusammenfassend kann man sagen, dass die Hauptelemente der Struktur der 
Spaltflåche von InSb-Einkristallen die Splatstufen und die Gleitlinien sind. Die Spalt- 
stufen, die in der Richtung der Spaltverbreiterung verlaufen, entstehen infolge der 
Schraubenversetzungen deren Verschiebungsvektor senkrecht zur Spaltfläche steht, 
bzw. in allgemeinem Fall die Spaltflåche schneidet. Die Gleitlinien kónnen unabhängig 
von Spaltstufen sowie auf den Spaltstufen beobachtet werden, wo sie gewóhnlich unter 
grossem Winkel (nahezu 90?) zum Spaltstufenverlauf stehen. Es sind zwei Gleitlinien- 
systeme, die von zwei verschiedenen Gleitflåchen herrühren, beobachtet worden. 

Infolge des Zusammenhangs zwischen Spaltstufen und Versetzungen sowie 
der Beobachtungsmöglichkeit von Gleitlinien auf der Spaltfläche, kann die Unter- 
suchung der Spaltfläche eine Methode zur Untersuchung der Defekte in Einkristallen, 
auch von manchen wichtigen Halbleiterstoffen darstellen. > 
Die Autoren danken Herrn Professor Dr. L. Sosnowski für die cuit Dis- 
kussion. 
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The Bloch spin wave theory of ferromagnetism is presented. Starting from the 
hamiltonian in the Bohm — Pines representation, the author obtains the Slater — Bloch 
equations and formulae for mean magnetization for three types of cubic lattices. All the 
calculations are made on assumption that the atoms are monovalent. At the end of the 
paper the value of the parameter A,, introduced in preceding papers, is estimated for 
crystalline hydrogen. 


1. Introduction 


In part II of this paper we have derived the hamiltonian in the BP representation, 
the appropriate subsidiary conditions and the dispersion frequency relations for electrons 
moving in a crystal medium. It thus becomes possible to apply the spin wave theory 
to describe the magnetic behaviour of cubic lattices. 

It is obvious that in our case the Bloch spin wave treatment must lead to results 
different from those obtained by Bloch. This results from the application of spin wave 
analysis to the hamiltonian, which includes terms with short — range Coulomb interac- 
tions. We must expect therefore a change in the value of the exchange integral and 
consequently a change in the value of the magnetization. 

All these considerations have been restricted to a model with one valence electron. 
Nevertheless the generalization to more than one valence electron does not offer any 
serious difficulties. As to the interactions, we have assumed that these exist only 
among nearest neighbours, since the particles act through the screened Coulomb 


potential. 


2. The Slater — Bloch equations 


. We rewrite the hamiltonian, the subsidiary conditions and the dispersion 
relations from part II putting Z = 1, which implies n = N, 
(295) 
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Ag Ne? li a , , 
+5 ho, (ala, + 1/2) T Nå m to= Fç. (2.1) 
[A «A, 2 [å] «A, 
N 2 
— > 00; > dn '« 
exp (i Å Tj) Io = 0, AKO, lA X A, (2.2) 
> op GÅ TN VU Az 
er a E ee Trey (2.3) 
= . mA? , 2.: 
w, = Op + 2 mio, 2 (A: p;)? + Bane |A| < 4e 
2 [sinc | 
F@) = |1—— | dr], (2.4) 
9 
m lin tA ' 
= wer ke ti A MA 
m* ( Fr): Gs) 


where m* denotes the effective mass of electron, T Pi» e are the position vector, 
momentum and charge of ith electron, respectively, Å is the Planck constant devided 
by 2z, V,(r;) = Vil, — ra) the potential energy of ith electron in the field of æ th 
lon — core, HE a, — creation and annihilation oscillator operators, w} — electron 


% 
An Ne? s] == (AnuNe*|mV )^. Langmuir fre- 


mV MV 
quency, A — reciprocal lattice vector, Å” — screening radius, V — basic cell 
volume, N — number of lattice points in the basic cell. 


When the subsidiary conditions are transformed to the form 


y exp (iA -) [L-E G B2mo, --hAt[2mosf + ...] p=0, 14 0/1 < 45 (2.6 
> Di 


dispersion frequency, oO, = | 


_ 


-. 
Il 
M 


it can be seen that (å "B;/mo, + hå?/mo)? is a small quantity if hA2/mo; <1, 
since, as we have assumed in part II, (A : D;|mo )2> Ay <1. We have used the 


former inequality when transforming the dispersion formulae. On the other hand 
N 


the expression > exp (i Å» F;) is always small (random phase approximation), hence 
i=1 


the subsidiary conditions are satisfactorily satisfied by arbitrary functions. All these 
reasons force us to choose the trial function p in a form which guarantees that the 
operator equation with hamiltonian (2.1) will be satisfied. The appropriate function 
9 for the ground state of the system can be chosen in the form of a product of ®, 
(ground state oscillator vector) and a function describing the behaviour of SRE 


Or isa 2 
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Before we write this last function in explicit form let us simplify the Hamilton opera- 
tor (2.1). For this purpose we observe that 


HI |A| <a, 
and 
p ho» (az, + 1/2) ©, = (1/2 på ho) ©, = É hoy + 
[A] <a, VEDR [A] <a, 
N fe 
3h ris. h? | hao, V 
ES b Add 
yi 4Nm?o, 2: > CP l6m?o, > Ag An? ur 
NA, [al <a, 5 


2 
N [4 
hV Ae 1 ney i SAR 22 
zh 8Nm?o,z? x Pi [ Er, 32m?o nz? AT Je a 
= 0 0 


h3V 


7 
+ ——— — = A] D. 
40 N Der Ta Vai T Amt, UL 4 ) ° 


Taking into account the transformation of these terms we obtain the changed 
hamiltonian in the form 


= Y pian + © Sg -F (he ry) + Då Von) Fr |F, — Fa) + 


150 


2Ne? ho V —3 RAV 
zo x he tas 3la a ^ T 224m, 


x, Hyp= Ev, (2.7a) 
or f w 


N 
ie = Y Bm 4 5 Y — Mica) + Y FRA), C) 


i=1 1255 ia LÉ 
# ç = Eg, ij (2.7c) 
, 2Ne? MNT BY. ur 
Era la e 


: AR ) (2.76) 


, V 3 9 
time = im (1— TRE 3-315 ^ TN Amon he 
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We can now write the part of p which depends on the electrons: 


pı (f4) NS,» CAT (TN) NS,» (on) 


ER 1 
AUDE" 


+ 


: (2.8) 
on (F1) NSN? (23). +++ on (Fn) sn: (On) 


where we designate by h; ..., h, the spatial distribution of r right spins and Ns, (7) stand 


for spin functions. We have (N) such functions  (h}, ..., h). The monoatomic wave 


functions p, (r;) are the eigenfunctions of the operator 


$99, = Ey Po X9—j)2m + V, (D. (2.9) 
As to the integrals 
eC i Pa (À) Pp (È) Pal?) på (1) dr dr, aA p. (2.10) 


we assume that they are small and therefore negligible. 
The resultant wave function for the whole system will have the form 


p=% X ah... h) øk... h), (2.11) 
(h^... hi) 


where the summation extends over all (N) distributions of right spins. The ah‘, Rs h) 
are the spin wave functions. We have to remark here that the description of the 
system by means of monoatomic functions is valid only when the screening is not 
too strong in the space of the elementary cell. 

Using the expression (2.8) we get 


Hp=% X a(h,..., h)Z e (hi, ..., h) 


(h^, D 
= a 3 A (2.12) 
= EO, >) gt sh) UR... h). 
(h^... Ri) 
or in matrix form 
å > a (hy, ..., h.) (hy, ..., h, EA h,, ..., h) 
(h^... hp) 
R- Ls cy qe + (2.13) 
se E > a (h,, ..., h,) (hy, ..., h |h,, ..., h,). 
TE | 


Following Bloch, we take into account only the terms which are connected with 
schemes of distributions as for example: 


x° rcu Fee < < 
PUOI TUE: a a+l N-1 N yh 
ac EM. pec al Mer (2.14) 


h visis h e e € ° ... ° ° 
LÉ OI LET wag @ at N-1 N 
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Hence we get immediately the nondiagonal matrix elements in the form 
(h,, x mh FADE? ies T (A), (2.15) 
where the exchange integral 7(A,) is given by 


Ha) =3 | v O 9 Open Ea FY 20409 ren] + 


+ Ben 0 Yen ©] + [26.0 re] + 


+ pom ka TE p e| + 27 (Ë —#p] dr : di”, (2.16) 
with 
Adv — 7, | 
2 in T 
Ga (r) = Va (z) Bl = ar] (2.17a) 
° 
alr- | 
e? 2 sin T 
PT NN m n = 3 2110 
y (T —* |) T$ ah - ir] (2.17b) 


The diagonal matrix element has the form 


(h... h. | |. Me uu N f ee V (r) di + 
V 


+E SË G, () d + 1⁄2 D [| 22 @) V (|F — |) g P) dø dø — (2.18) 
a BV a+b V 


(MND | pari © pe © pari (P) pa ©) | Ë Ca) — a V, an 
y 


+ Ë Cats 0) — 37 Ver o + |: G, (P) — Ve e| L Bes (n) = 


— Var e»| 2. arin — ha di? = — NE + 6(9 — (N, + Ny) J 09. 


We have designated by @(A) the Coulomb integral 
ewan fr SIE (0— 2 a+ y fo (P) Gp () di + 
a+b y 


txux[feogrvü- Pos @) ar or _ (2.19) 
“G= 8 V 


\ 


and N,, N, are numbers of pairs of right and left spins. 
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Since we have assumed that the overlapping integrals are small we have 


Ev > 
(hos. h, | hi, ..., hi) = ài... CA ME 


which together with (2.15) and ra leads eventually to the following Slater — Bloch 
equations 


adv) (h | £ | h) = [E— (h| £ | B] a (h), (2.202) 
hh) 
or 
JA) D la (h) — a (h) + e (A) a (h) = (2.20b) 
h' GER) 
where 


e (4) = E— (h || B) + g (A) = E— — NE, — (4) + 
EM EN 1) 702. (2.20c) 


E 


We have designated by h; h’ the distributions h,, I h,; hi, wh 
Since, according to the proof of Dyson (1956) the spin complexes of Bethe (1931) 
do not exist in low temperatures, the general solution of equation (2.20b) will be 


a (hy e h) =C > (00) I exp (io - h), (2.21) 
Po) s= 


and the energy of spin waves 


e) = 709 X XL exp (a9 4), (2.22) 
Ul 


with the restriction r < N^, Here C is a normalizable constant, A the vector to the 


nearest neighbour, the sum over À indicating summation over all nearest neighbours, 
^ 3 er es 
P(w%) is commutation operator, w® = 2zJN% >} s® b°; b? are vectors of reciprocal 


p=1 
lattice, Ou = 0, + 1,..., #(N% 9/2 


3. The spontaneous magnetization of the system 


Introducing a magnetic externat field of strength H directed along the positive 
z — axis we get in the hamiltonian the supplementary term 


N 
>, + Eh Fe >, >, 
(Ri, «5 Ry | Ne ôi H |hi, ..., hv) = Hu [r — (N—n)], (3.1) 
= 
where we have taken into account the fact that 
oi Ns. (o,) = + CARTE ng (o;) for right spins, (3.2a) 
oi ne (0) = — 9,59 m; (0; for left spins; (3.2b) . 


u is Bohr magneton. 
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Keeping in mind that spin waves obey the Bose—Einstein statistics we obtain 
from (2.22), by (3.1). 


£o. ..., e Fear aoa, (4) = — NuH + 


+ 2) ng [eg (4) + 2H} (3.3) 
P 


3 
E TAG cm 2 > 
where 3 ng = r, eg (A) = Q0) V [1 — exp (667 ` A), 6; = Fy y, aÀ 


? å i=l 
and p; are numbers from the sequence 0, + 1, ..., 4: (N^ — 1)/2, but taken in such 
a way as to avoid repetition in EX (Ay), un = V 2,2... 


The partition function will be 


N 
Z = exp (NHu|kT) V. Y I] exp E T (eg (A) + 25g) | GER 


r=0 (n$) P 


hence the average magnetization 


== 9 InZ 1 
M (H, T) = kT NIE 
( ) oH A DD |z | (3.5) 
SHI. 


$^ exp 


The series på exp |- IT [e AES 2n): in (3.4) is convergent if ez (À) > 0. In 


"p 
connection with this inequality we get J(A,) > 0 as a condition of the ferromagnetic 
behaviour of the system. 
Since the calculations are quite identical with those of Bloch we give here only 
the final results for the spontaneous magnetization of three types of cubic lattices: 
a) simple cubic lattice 


M, (0, T) = Nu p- (25) '| r= Di wane (9-409 69 


b) body — centred cubic lattice 


ær 3/3 
TA (0,7) = Nu |! . [s ss) ] 99-102 AG G 


c) face — centred cubic lattice 


5 IR f 
M, (0, T) = Nu er im å = = i Od) = J 9l. | (3.8) 


The exchange integrals are computed for the nearest neighbours, with the excep- 
tion of the body — centred cubic lattice; 7,(A ) denotes the exchange integral between 
the nearest neighbours, and /7,(A) between the next nearest lattice points. 
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4. Application of the theory to metallic hydrogen 


In order to compare the Bloch analysis and the theory presented here, we apply 
the results obtained to the hypothetical metallic hydrogen. In this case the wave 
function, the ion — core potential and the energy of the electron are given by 


1 % T c? e? A2 
T) = (— —— |, Vn=—-,E = — 5—, % = — 4.1 
o@) = (a) vl.) mine: Gl) 
Since we are using the monoatomic functions instead of the Bloch ones, we must in 
à T V 
all integrals pass to the limit V — oo, N — co, — < oo, and compute the energy 


for one electron. We have from (2.19) 


1 ip [í 2 has ER ' 
V 


Agr 


2 [vol a) roai foh- 
V 0 


a£B y 


edt (4.2) 
2 sin T Nin. MM 1 
SES iF = ir] V (r) dr +57 IK zes 
0 a£B y 


alr=r| 
2 sin T vp MON 
zc J = ar] ps (^) dr - dr”, 
ó 
and the exchange integral, according to (2.16), will have the form 
Alr—ral 
4 Fr P^ an P. 
F(A) = Ja—— f Parr Å) Pa (P) Pari) Pa (7) |". O) ih EE dede 
y o 
Aglt—re 411 A. 
e 
Va a RERUMS ZW ug itas dø -di’ + (43a) 
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TEKTI hA Å E z > 
tay E - ET [ erc eo + 
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where we have designated by 
(4.3b) 


OB = fous (T) Pa (F) Yarı (È) Pa (1?) |. (r) + Vous (r) + rum ix dr : dr’, 
V 


: . : e 
the Bloch exchange integral without the ion — core Coulomb energy —, xP. 
Tap 


Let us suppose that all the spins are directed along the external magnetic field. 
Then we have from (2.20c) for the energy per electron 


E Nc 282 LME € 1 s h V 2 k 
Pes + pay (Sto, +2) Le DOSE MUT 
h? % 7 T 5 E? 
ARE A AS (4.4) 


The value of the screening parameter À, must be obtained from the condition of extre- 
mum of energy. Hence 


% [E 2e? Im e2\ 22 i Å 
dA, a CU rms rt Jr ILE ` 2: ah 
hs : 
en! Ua 6 09 — 57 JU =0, (4.5) 


where S denotes the number of nearest neighbours in the basic cell V. When we 
replace Å, by À, = À, * ag and å by ô, 43 = ó, 473 a$, where A’ = zs À stands 


0 
for the distance between two nearest neighbours and ô, = 1, 1/2, 1/4 for the simple, 
body — centred and face — centred cubic lattice, respectively, we get 


, 2 23 2 6, 2 
xls) Rap Ed soie eure Qm RU LEUTE. 


9 AN G Oa AU ms TN M 
ko Ag 
= Ae) = Q). 4.6 
32m? wp ag x? "t PL €) = + 5 IA) = (46) 
1 9) ME. š 
The exact calculation of — 5 P where only the nearest dan. are 
taken into acconrt gives 
Lu e? ae qe hae: ~ has 2 a 
N 31, £9 no ID + Ig Amw, a2 g 


9 
For the expression 277 (A) we get - 
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gi JU = FE ep (24) rar a) Tr 
2 cuim] | er E EE — 
nn + | i | + Gp 
Aurea) =. Me tau 


We shall look for a solution of the equation (4.6) under the assumption that 
2 

<A, «£l. This condition agrees Mer Wi which we have already used. 
44 mo, 
Figure 1 shows the dependence of Å, on A’. 


| : f tag 9 [E 
In Table I are given the numerical values of the derivative —- — | — 


e JA \N 
for A Li 
Table I 
simple cubic lattice body-centred cub. lattice face-centred cub. lattice 
A Tao 9 E 
~ e294, NN 
RE SN à 
2,03414 0 0 
2,035 +0,0396—0,2309 22 | 0,4141 | 4+0,0514—4,4619 22 |0,1073 | --0,0770— 6,8563 22 |0,1060 
2,040 4-0,2702—4,6871 42 |0,2401 | --0,4054— 7,2028 42 |0,2372 
2,045 +0,5006—4,9185 42 |0,3190 | +0,7510—7,5578 42 10,3152 
2,050 +0,7200—5,1552 A2 |0,3737 | +1,0800—7,9158 42 |0,3614 


0.50 


0.25 


2.085 2.040 2.045 2.050 4’ 


os] 


Fig 1. The s-/ curve represents the dependence of À, upon relative distance A’ for simple cubic lattice, 
b-l the same dependence for body-centred and J-L for face-centred cubic lattice. 
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Ås can be seen from Table I, the assumption about all spins are directed along 
the external magnetic field (saturation magnetization) leads to a maximum of energy 
for the system (the second derivative is negative), while only the minimum value of 
energy would make the system stable. Hence we conclude, that metallic hydrogen 
can not be ferromagnetic. | 

Since we have already obtained the screening parameter À, = À,a, for a range 
of lattice interdistances, we can now calculate the value of the exchange integral and 
compare it with that of Bloch. Integrating the expression (4.8) for a given À, we get 


2c? P P il , 3 8— 3 T , 
dh) = a. tee a I tian ær | 442 4A Tr 
Teige ,, 32+ A) 1 1 
m sin Aa A «| 258971: Ag lm sips) | ene d 
Ag” 
6(1+ 2A, 3 Å x sin T 
4a E VE ^m + VE (1 -H 2A d- Å 2) ies ihe =—= (4.9) 
Beaded | dcc dU aet a] + 
3! 37 Int 
T. hô, (1 2 A^ ) "m 
20zt mo, ao 


On the other hand we have, according to (4.3b), 


8 LEO e m. 49 h s > TTE 
a E me are p (C + hn) 
6 m QT. 
x sts tung des exp (—24) + |z — 4 + 
| Hot ns 
OU nå E; (— 44°) exp (24’)/ 4^ — d d 
ae = ar) bac 2014), C = 0.57721566, (4.10) 
where E; (x) is the integral exponent 
— E, (— x = ET x > 0. 
T 


The comparison of our (A) with that of Bloch gives the Figure 2: 

We see from Figure 2 that the values of our exchange integral do not exceed 
the Bloch ones, therefore these are always negative. We have thereby obtained once 
more a proof that metallic hydrogen cannot be ferromagnetic. 
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Ac) - 4:9 
A (Ae) = EG in% 


2.030 2.035 2.040 2.045 2.050 


Fig. 2. The mean deviation of the values of the exchange integral from the similar Bloch quantity: 
As—ı (À) for simple cubic lattice, A,_7(A.) for body — centred cubic lattice, Ar (À) for face — 
centred cubic lattice. 


I consider it my pleasant duty to thank to Professor S. Szczeniowski for suggesting 
the subject of this paper and for many valuable and helpful discussions. 
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THE GENERAL FOKKER ACTION PRINCIPLE AND ITS APPLI- 
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À method of constructing a Fokker action principle (a principle in which only the 
dynamical variables appear) for the equations of motion of dynamical systems is given. 
This method is used to finad the action principle leading to the post-Newtonian equations 
of motion of a perfect fluid in general relativity theory. The Lagrange function appearing 
in this principle is discussed for the case of the motion of a fluid concentrated into drops. 
The group of isolated drops of the perfect fluid was introduced as a miodel describing the 
motion of the system of bodies in general relativity theory. The behaviour of the lagrangian 
as the drop dimensions tend to zero is investigated and the renormalization of the equations 
of motion of point singularities in general relativity theory is discussed. The calculations 
for finding the equations of motion in general relativity theory by this method are much 
simpler than those by previously known methods. 


The general theory of relativity is thus far the only theory in which the field 
equations embrace the equations of motion. The equations of motion for bodies of 
large mass were introduced into this theory from the field equations for the first time 
by Einstein, Infeld, and Hoffman (1938), and subsequently by means of other methods 
by Fock (1939), Einstein and Infeld (1940, 1949), Papapetrou (1951), and Infeld 
(1954, 1957). In addition, a variational principle leading to the post-Newtonian equa- 
tions of motion for bodies was given and discussed in the papers of Fichtenholz (1950) 
and Infeld (1957). 

The motion of bodies has been discussed thus p on the plane of the field equa- 
tions and equations of motion, but the possible connection of this problem to a field 
action principle in general relativity theory has not been discussed. The approach to the 
problem of the motion on the level of a variational principle has allowed us to give a new 
method for obtaining the equations of motion in general relativity theory for a rather 
wide class of dynamical systems satisfying the conditions given in the next section. From 
the point of view of the calculations the method is more economical than the methods 
used till now. The method discussed here, however, cannot be used in considering 
the relation between the field equations and the equations of motion. If, therefore, 


(507) 
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one is seeking an effective form of the equations of motion of some dynamical system, 
then this form is most easily obtained by the method given here, in which the problem 
is treated on the level of a variational principle. If, however, one desires to consider 
the relation between the field equations and the equations of motion, then this can 
be done only by the traditional methods. 

In the method presented here, one constructs the Fokker lagrangian (i.e. the 
lagrangian in which only the dynamical variables appear) of the dynamical system 
under consideration whose Euler-Lagrange equations will be the equations of motion 
being sought. The procedure is thus somewhat different than e.g. in the paper of 
Fichtenholz, where in the case of a special dynamical system — a system of material 
points — the Fokker lagrangian was obtained by trial from the form of the equations 
of motion, while in our case, the Fokker lagrangian is constructed from a field action 
principle without knowledge of the equations of motion. This procedure can be used 
for different dynamical systems, and not only for material points. In the case of the 
motion of material points, there exists another method of constructing the lagrangian, 
which was given by Infeld (1957); this is not, however, a consistent field method, and 
in its application encounters difficulties which must be by-passed in a rather artificial 
way by introducing an alternate field g% g. On the other hand, the calculations connected 
with this method are more or less of the same degree of difficulty as, e.g., the calcula- 
tions by Infeld's method (1954). 

The simplicity of the method presented here results primarily from the fact 
that, at present, in order to obtain the equations of motion it is sufficient to know 
the field in a lower order of the approximation than was the case previously. For 
example, in 1938 to find the post-Newtonian equations of motion it was necessary 


to know the quantities oo Aom hoo, Rom and Apn; in 1954, in the paper of Infeld, it 
RT. ge maga: 4 


2 
was already sufficient to know oo, how Nog; at the present time it is sufficient to know 
et 
the simplest fields hog and Aom. The effectiveness of this method is illustrated by the 
2 3 


fact with its use, after relatively simple calculations, the post-post-Newtonian equations 
of motion of material points have been found (Ryteñ, in preparation). To give the 


equations in such a high approximation it sufficed to know the fields hoy, hom» hoo, Ryn 
2: Fg sand TE 


while it would be practically impossible to find this approximation by the traditional 
methods. 


1. The Fokker variational principle for a system of material points. 


In this section we shall consider the field described by the system of functions 
y™ (x^) (M and the next capitals in the alphabet run over the set of indices of 
the tensors, spinors, etc., u = 0, 1, 2, 3). We assure that this field interacts with 


` ` a ; 
points sources having world lines £* (x0) (4,B,C = 1,2, ..., p. where p is the 
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number of particles, while for the parameter of the world lines of the particles we 
A 
0. de" 

take simply the time x®; also, we assume that B "ri = 1 and we assume that 


the dynamical laws of the above system of particles and field result from the 
principle ôW = 0, where 


A X A A x 5 
WE, pe] = f do? A (E, på) + f d? [age LOM, g); and) (1.1) 


A A4 
y" = f dax ô (z = Ë) y (z, 28); 
A SANA 


y= Exe vt ree SETENE SS), 


The generalization of the considerations of this section to the case in which the lagrang- 
ians À and L depend on higher order derivatives does not present any difficulty. 


A 
We calculate the variation ôW produced by the variation of the world line 6&4(x®) 
and the variation of the field óyM (x^) 


re far de fao fan, 
2: I D (a9) RA i ST 


ace 


(1.2) 
Ate A x< ° 


SF PE da 


L+ J d,o Ph öyM Ny» 


by Dj; we denote here the hypersurface ee the four-dimensional region 
contained between the hypersurfaces x? = x'? = constant and x? = x''? = constant; 
n, is the four-vector normal to $7; also, 


u 
óW oA 4 d 8A = 9A 
A a a oa Dat a (1.3) 
ge Qo SER DE 984, 
OW v Spy Lay n r VN 
dy (v^) T By Da BM, Sy, 


4-1 09M 


'The indices characterizing the transformation properties of the quantities considered 
(e.g. the indices M, œ and a) are subject to the summation convention, while the indices 
numbering the particles are not subject to this convention. 

If, in accordance with the action principle, we assume that the variations of the 
field dy" vanish on the hypersurface > and that the variations of the world lines 


A 3 A 
1) The function Ô(x — 8) = = IT ó(xi — &) does not denote here the ordinary Dirac function, 
i=1 ` 


but the modified ó fonction introduced by Infeld and Plebanski (1956, 1957). 
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of the particles vanish at the ends of the interval < x"°, x’ >, then from the principle 
ôW = 0 we obtain, by (1.2), the “field equations": 


A 
OW [£^, y^] 


OU 1.4 
ye T be 
and the “equations of motion“: 
B 
LAs DE HMM (1.5) 
029 (29) 


Using the terminology introduced by Infeld and Plebañski, we call equations (1.5) in 
which there appear the unknown fields y the equations of motion of the second 
kind. Let us now assume that the field equations (1.4) can be solved for an arbitrary 


A 
motion £^(x9), and let us consider one of the solutions of the field equations 
A 
y = y™ (x^) [£*]; (1.6) 
such a solution will, in general, hats functional of the world lines of the particles 
A 
&°(x0). The solution (1.6) obviously satisfies the field equations identically, i.e 
A B 
oT [£*, y, [£*]] 
dy” (x°) 


If in place of the fields appearing in Eq. (1.5) we substitute the solution (1.6), then 


= 0. (1.7) 


A 

we obtain a system of equations in which there occur only the variables £^ (x°) and 
their “time“ derivatives. In accordance with the terminology on Infeld and Plebañski 
these are the equations of motion of the third kind?). In symbolic form they are 

A B 
OW [š#, på [E] 

A 
à (84 (4^) 


where we have given the argument of the functional derivative in the parenthesis 
to stress that we have to do here with functional differentiation with respect to the 


= 0; (1.8) 


?) According to this terminology the equations of motion of the first kind are the equations of 
motion of a test particle; the equations of motion of the second kind are the equations of the form (1.5), 
ie. the equations in which there appear the unknown fields yM (x^); the equations of the . kind 


together with the field equations create a system of equations with the unknowns y™ (x4), Eat): the 
equations of motion of the third kind are the equations of the form (1.8), i.e. the equations written only 
with the help of variables referring to the world lines of the particle. a equations, for suitably formu- 


lated initial conditions, allow the determination of the world lines fa (x9). 


1 
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A A 
argument £^ (x°) occurring in W[E*, y”] explicitly, and not through the intermediary 


A 
of på (x^) [°]. 
We shall now explain what we mean by the action principle of the Fokker type. 


A 
Such a principle is a variational principle ô W,=0 where W, = Wp [E°], from 
which the equations of motion of the third kind (1.8) result directly. There are physical 
theories for which such a principle can be formulated and those in which such a prin- 
ciple does not exist. For example, for the Newtonian theory of gravitation or for the 
post-Newtonian equations of motion in general relativity theory (see Infeld 1957) . 
such a principle exists. 

Thus far, where the Fokker et existed the action was obtained by trial 
from the specific form of the equations of motion of the third kind. We now want 
to consider whether we can find a more general criterion which would allow us to 
determine whether there exists a Fokker principle in a given theory and whether 


A 
one can then find W,[£°] without knowing the specific form of the equations of 
motion. 
Some suggestions in this direction can be obtained from the shape of Eqs. (1.7) 
and (1.8). If in the world of functionals there were to exist an analogy to the theorem 
of the differentiation of compound functions 


d _ of , Of dy 


then we would have 


B G 
ôW [é pr [EY] — SW [&, y REI ; TAG y" [£z] D SIO 
esu cola nsu SP Gr) Fidel de 
DES (a9) Ô (e (x°)) 6&4 (x?) 

and therefore, by (1.7), 

B G B c 

OW [&, vr EJ] _ óW E, vr ET, 

A ber A 

de (x?) à (84 (x*)) 
Eq. (1.9), in general, is obviously false, nevertheless the above considerations suggest 
that if the Fokker principle exists, then one should seek the action Wp in the form 

A B C 
We [E] = W |ë, y", ET, (1.10) 


obtained by eliminating the field from the total action (1.1) of the given theory after 
| inserting into it one of the possible solutions of the field equations. 

One should, of course, check which conditions must be satisfied by the given 
physical theory in order that the hypothesis expressed in (1.10) be correct. To do 
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this we calculate the variation of the functional Wp given by (1.10), 


å » À 
TE J w Y "i. yM[E]] je 1 Jefi ó] [E y OW [Ë, pr [E] TM 
f£ ee») dyM (x) 


gx"? 


+ | dro Ppa yna 


Ë A A Vh 
+ Y Da 6&4 
At x 


A 
where ó*yF are the variations of the functional y? (x")[£^] resulting from the 


+ [do Pre d*yR ng, (1.11) 


o=x'T 


A 
variations of the world line £^ (x9), that is, 


syr Ÿ ar EN så oa, 
A=1 ye öga (x^) 


A 
while p; and P/* are expressions which can be obtained from (1.3) by inserting into 


A B 
the formulae given there for p, and P the field y (x) [£^], o is a hypersurface describ- 
ed by the equation x? = constant, while v is a three-dimensional time-like hyper- 
surface (at infinity with respect to the variable x). 
The second expression in (1.11) vanishes because of the field equations (1.7). 
We Rosse see that the action principle dW, = 0 for arbitrary variations of the 


world line "m vanishing at the ends of the interval < x”, x” Es leads to the equations 
q 


of motion of the third kind (1.8) if and only if the field y? (x^) [£] satisfies the condition 


a 


dio. P; ORUN n, 2: dør På ó*yf n, = 0 142 
pn y” fuis å (1.12) 


A 
for arbitrary variations ó£“ vanishing for x9 = x'9 and « = x’, 


A 
Condition (1.12) is greatly simplified if the solution z (x^) [£^] is not a functional, 
but a function of the world A that is, if 


v" (^ G = y (x, Ë S» & Mu D (1.13) 


Then ö*yR is a linear function of e variations D (x), A (49), ... and the require- 


ment that the variations DEA), E (x), ... vanish for x° = xx" implies the 
vanishing of the first expression in (1.12). On the other hand, the second integral 
in (1.12) can be written in the form 


x"? Ñ 
Sedo Pi 0*y m, = f dx $ do Pre ó*F n,, 
+ x EL E, 


where F is a closed ordinary two-dimensial surface at infinity. Thus, for solutions of 


BR ii ce 
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the type (1.13) condition (1.12) passes over into the requirement that for every a? the 
integral 


Pda PF open, —0, (1.14) 

A A 
for arbitrary variations ó£^ (x0), O (x9), -, vanishing for x9 = x’, x’. Condition 
(1.14) is e.g. satisfied for those fields 4? whose asymptotic behaviour at infinity guar- 


A A 
antees that for arbitrary variations ó£^ (a9), am EE 


På ô*yP SG) 
at least as rapidly as 773; r = (a? 29)%, 

Within the framework of the considerations presented here it is difficult to say 
something more about conditions (1.12) or (1.14). It should be noted here that in 
a given physical theory there always exists an arbitrariness of the choice of diffe- 
rent solutions of the field equations. These solutions describe fields of different 
physical properties and have different asymptotic behaviour at infinity. (As an example 
one may cite the solutions with retarded or symmetrical potentials in electrodynamics). 
For some solutions condition (1.12) may be satisfied and a Fokker principle will then 
exist, while for others it will not be satisfied. Physically speaking we should expect 
that the Fokker principle will not exist, and therefore condition (1.12) will not be 
fulfilled, for "radiative" solutions, where there appear in the equations of motion 

A 


damping terms containing odd derivatives of £“ (x°) with respect to x°. On the other 


A 
hand, in the case of "standing waves“, where there appear even derivatives of £^ (x) 
with respect to x? in the equations of motion of the thirs kind, or, more accurately, 
when the equations of motion are invariant under the inversion x? — —x®, one should 
rather expect the Fokker principle to exist, and therefore conditions (1.12) to be 
fulfilled. 


2. Fokker action principle for a dynamical system in general relativity and the EIH 
approximation method 


We shall now consider the formulation of the previously mentioned idea of con- 
structing an action of the Fokker type to find such an action for the equations of 
motion of a dynamical system interacting with the gravitational field in general relativity. 
In the following sections we shall concentrate chiefly on a special dynamical system, 
ie. a relativistic perfect fluid?) For the time being in the general considerations of 
this section we shall not single out any dynamical system. Thus, for the present, the 
variables y? (x) will be the components of the metric tensor g,,, (2%), while the variables 
describing the dynamical system under consideration will be denoted by dy 


3) Other dynamical systems, such as systems of point particles represented by singularities of the 
field and a system of point charges, have already been discussed previously by one of the authors (Ba- 
zanski, in press). Å system of rotating bodies has been discussed by R. Michalska (in preparation). 
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(M = 1,2, ... f, where f is the number of variables). In the case of material points 
A 


the variables dy will be the coordinates of the world lines of these points £^ (x?) and 
their derivatives &% (x°), while for liquid the variables dy; can be the Euler velocity 
field, pressure, and density. In the previous section in order not to complicate the 
solutions we employed a system of material points, since for different dynamical 
systems the action principle is formulated in a different way. 

We assume the following action in general relativity theory: 


W= = ji G (Bun Euro) V-g dax + j A (dm, Sun Surio) V-g dax, (2-1) 
where 
G = g% (dø) Gu) — lan) Usb (2.2) 
g = det (gag), 


and A (dm Sep Zapio) is the Lagrange function characterizing the dynamical system. 
From the action principle with respect to the field g,, we obtain the Einstein equations 


R» — = gi? R = — 8x T”, 


re et 2 UM, (2.3) 
y—a 9g, 9x" ILuvio 
and from a suitably formulated action principle with respect to the dynamical variables 
we would obtain the equations of motion of the second kind for the dynamical 
system. If we would wish to make use of the considerations of the preceding paragraph 
and construct the Fokker action Wp[dy] for a class of solutions of the field equations 
(2.3) satisfying condition (1.12) and then find the equations of motion of the third 
kind, then we would encounter a basic difficulty. This is due to the fact that we cannot 
solve Eqs. (2.3) for arbitrary motion dy; since the condition of the integrability of 
Eqs. (2.3) is the satisfying of the equations of motion by dy. This difficulty can be 
overcome by use of the EIH approximation method. We shall recall here the basic 
assumptions of this method. 

We assume that the dynamical system under consideration is an isolated system, 
i. e. the curvature tensor of space-time vanishes at infinity, and we assume there a Gali- 
lean reference frame: 


lim Sap (x6) = Naps 
|x| +00 (2.4) 
Noo = 1, Nom = 0, EN F Os’ 
for large x^ (a = 1, 2, 3), i.e. when the space of variables x^ approaches a Euclidean 


` 1 . w .` * 
space, the "distance" | x | = (x* x‘) e passes over into a Euclidean distance, while 
x"/c passes over into Newtonian time. 
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We next define the quantity h,, (not a tensor) 


hag Bat lap» (2.5) 


which is not riecessarily small; it must, however, in accordance with (2.4) satisfy the 
condition 


lim hy, = 0. (2.6) 


|£| ^ oo 


We assume further that h,, is an analytic function of the parameter c7! (cis the 
velocity of light at infinity) and we represent h,, in the form of a power series 


hap = hap a; hap Åp hap Tee Cee 


where te, = hr, ch and h“) does not depend on c. We assume that the derivatives 
i 


2 
io^ NE 
Oxk P 

I : Å 
239 hg = Papio are proportional to c^! 
i i 


the index zero, e.g. hago: The physical reasons for this treatment have been given in 
Eg 


hapik are (as is the case for h,,) proportional to c and the derivatives 
i i 


, which we denote formally by writing 1 under 


several papers on the equations of motion. They will also be found in the monograph 


by Infeld and Plebanski, "Relativity and Motion“. It can be shown (see Infeld and 

Plebanski "Relativity and Motion“) that from the condition of correspondence for 

c — oo between the equations of motion of a test particle in general relativity theory 

and in Newtonian theory it follows that hoo = Aom = 0 in the expansion (2.7). According 
1 1 


to the general principle of relativity the quantities h,, are not, of course, uniquely 
determined by the field equations, but to an accuracy of an arbitrary transformation 
of the coordinate system x“ = x'* (x"). In the EIH approximation method, according 
to the assumptions about formulated above, one considers a class of transformation 
which: 


1) preserves the condition lim zg, = ap and at infinity passes over into an 
|2| +00 
identity, 
Ru 1 
2) does not change the type of expansion of the quantity p in terms of us 
Introducing the generator a“ (x*) of the transformation a'* = x'* (x^): 
cc oa) 


we can formulate conditions 1 and 2 in a formal way 


1) lim g (x) = 0, 
[2| + co 
2”) a = a? + a? + a + ..., 
TOEKS 4 


j ab esas qug sats. 
1 2 3 
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where a^ is proportional to c. Transformation of the reference frame generated 


by & fl % (a, = Nag a) produces a change in hg, starting with the orders Poe 
I+1 
Der pres ei 
HO HE ; > ’ 
h'o (x'*) = hoo (27) — 24010 (z), 
1-2 1+2 i+11 
hive (= ‚Ro n OY — aon (x*) — a, o (x^), (2.8) 
jie TES 1 1 


Re vin (20%) — amn (2%) — a, (x), 
1 l l l 


while the further terms of the expansion of h,, depend on a, , a, in a more com- 


I+1 d 
plex way. 
The group of quantities A,» ho, > hg therefore transform together under the 
l 1+1 I+1 
influence of the transformation generated by ao, a,. We indicate this by writing the 
IF1 1 


expansion (2.7) in the following way: 


EE Eg EN de lca aioe 
2 3 4 6 7 


Ya 
ho, = — 0 = hon RE hon "E hon sh hou ^ hon el hon pisa (2.97) 
/ på pa Å / 
adiit P A Á 
7 a 7 Ae Z 
Porn = Ryan cs Run d en gk hinn er xls Ryan T hn "T B 
1 2 3 4 5 6 7 


the dotted lines join the group of quantities transforming in the same way. In order 


to preserve for c — © the correspondence with Newtonian theory we put hoo = hom 
1 1 
= 0. Å more accurate analysis of the transformation properties of quantities h,, 


(details will be found in the monograph of Infeld and Plebanski) made for rather 
general assumptions about the dependence of tensor 7 “B on Zag Shows that by a suitable 


choice of generators dp, % and dp, a, (i.e. by narrowing the class of reference frames) 
SM 4 3 


the groups A,» Ron Roos ud Bn how hoo can be eliminated; as a result, the expant 
I5 dn E T 


ion (2.9’) passes over into 


_ 6 7 
ho, = 0 +0 + ho, +0 + ho, + hon + Pon +..., (2.9) 
5 6 7 
h, = 0 + hu +O +h, + hy, Fn + An + > 
: 2 - 4 5 6 7 


while the remaining groups of h,, written above cannot, in general, be eliminated by 
i 


The General Fokker Action Principle ő17 


a choice of the reference frame. In the expansion (2.9) there still remains the free 


choice Tolqa de... ap Ay As ... (in Appendix I we narrow down still further 
3 5 6 2145 


the reference frame fixing generator a,). We can, of course, always return to the form 
2 
(2.9) from (2.9) by means of the transformations generated by Ay = ao + ay, and 
25904 
a, = a, + a,. We note that in papers on the equations of motion one employs the 
1 


2 
approximation procedure of "skipping every second order“, that is, one assumes that 


hag = xs at for af = 00 and mn; 
hon = 2 hop 
=1 2;+1 


mn 


this kind of expansion is the same (2.9) up to the order Aq, hop» Amn» which allows 
EG 


the determination of the post-post-Newtonian equations of motion. 
We should also assume that the dynamical quantities can be expanded into 
a power series in c! 
dy = dy + dy +.» 
TM TM +1 
here ry is the lowest order term of the power series. If e.g. we take the velocity compo- 
A 


nent of the material point oi as se quantities, then 


A 
pe >> 4 + Ei d 


and in this case ry = 1. 
In the EIH approximation method only the quantities h,, are expanded into 
a series. The expansions in 5,5 splits up the field equations into the respective orders 


1 
and the individual h,, are determined from the field equations of the corresponding 


Li 
order. On the other hand, we do not expand the dynamical quantities into series, 
but formally treat them as quantities of the order ry, and determine them from the 
equations of motion which are correct up the order J, i.e. from equations which ensure 
that their solutions will agree with the exact dy if one neglects quantities vanishing 
at least as «#1. This procedure considerably simplifies the calculations. 

After this digression on the approximation method we return to our basic problem 
of whether the action (2.1) after elimination of the field g,, — which is a solution of 
equations (2.3) — passes over into a true action of the Fokker type for the dynamical 
system under consideration. Obviously, this problem can only be solved by approxima- 
tion methods. For, the approximation methods allow, as we shall see in a moment, 
the separation of the field and motion, and the individual approximations of the field g,, 
are, in general, a function of the dynamical variables and not a functional; we can there- 
fore make use of the simplified conditions (1.14). Our procedure can be outlined as 
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follows. By expanding the field g,, appearing in (2.1) into the series (2.9), we obtain 
an expansion of the action 


V= Ma IE EN 
1 5 


3 


We write the partial sums of this series 


på teu 
>3 1 3 
W=-WıW+W, 
+5 i: 3 5 


e e e o @ e o9 9 @ e alle n9 9 | nw | n 


(in the specific cases, in general W = constant) and in turn with an accuracy to the 
1 


order 3,5, etc. we should verify whether conditions (1.14) are satisfied; these conditions 
take the form 


& O* gap ny dao = 0; (2.10) 


we take for the surface of integration K (| x |) simply the sphere x, x, = constant. 
GV—g is obviously the homogenous quadratic form with respect to g,5, which we 


can write 
G y = = Cafuves EET PR (2.11) 


here C*#e? is the tensor density constructed from tensor Eag Components not con- 
taining the derivatives g,,),. Thus in order to check condition (2.10) we must 
investigate the behaviour at infinity of the expression 


I y—&C O* gap = 2CaBurck p. s O¥ gu, (2.12) 
9,5 Ë 
expanded up to the order being investigated. 

In order to eliminate from (1.1) the field g,, one has to avoid a fundamental 
difficulty connected with the fact that the conditions for the solvability of the field 
equations (2.3) are the equations of motion of the dynamical system which we wish 
to find precisely from the Fokker action. This difficulty is removed by the EIH appro- 
ximations. 

In order to calculate the action W we need to know hog and hyn = Smn Roo (see 

2 


3 2 2 
Appendix I); we find these quantities from the field equations for an arbitrary motion; 


we verify condition (2.10) with an accuracy to the 4th order and after eliminating from 
W the fields hoo and h,,,, we obtain 
3 2 2 


Wp [dy] = W la (dy), d]. 


This functional is the Fokker action from which we obtain the Newtonian equations 


mc 
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of motion permitting the determination of d,, to the fourth order, i.e. dy: From the 
+4 
field equations we can now determine ho,,, a function of the Newtonian4) motion du 
3 +4 
(see Infeld 1957 and Appendix I). We next check whether in the sixth approximation 


the condition (2.10) is fulfilled. Knowing hom (dm). we can write the action W, de- 
Pd 
pending on the Newtonian motion ath V, Ta, 
5 +4 


Wr "m W, (dy) + We [dm]. 
+5 +3 5 24 


In Wp, however, we can replace the argument dy by dy, understanding now that 
5 24 
dy is a motion which is arbitrary begining first from the post-Newtonian appro- 
ximation order; we then obtain 


-5 


Wr [du] = VEG [dy] rie ds (dul- 


The functional Wp [dm] is a good Fokker action (if in the sixth order the condi- 
+5 


tion (2.10) is fulfilled), from which we obtain ihe post-Newtonian equations which 
specify dy. 
-6 


In the present paper we have stopped with this approximation, since the above 
corrections to the motion are too small to ascribe to them any physical meaning. It is 
possible, however, to continue this procedure as long as conditions (2.10) are fulfilled. 
From the form of expansions (2.9) and from the analogy to electrodynamics one should 
expect that these conditions will still be satisfied for G, for this is the highest order 

8 


in which gravitational radiation does not appear. As we already mentioned, condition 
(2.10) for G was verified by J. Ryte for the case of material points described by 
8 


singularities of the field. Ryte proved, in this way, that there still exists a Fokker 
action Wp [dy]; Ryteñ obtained this action and the post-post-Newtonian equations 
-7 A 


of motion for £^(x9) with relatively little calculations (considerably simpler than the 
calculations leading to the post-Newtonian equations in the paper of Einstein, Infeld, 
and Hoffman in 1938). This fact proves that the method given here for finding the 
equations of motion is substantially more effective than the methods used previously. 


3. The action principle and the equations of motion for a perfect relativistic fluid. 


We postulate the action for a. perfect fluid interacting with a gravitational. field 


? 
ye | ae Ee le EN — J (8.1) 
ò 


4) Or more rigorously; a function of a motion which will agree with the Newtonian motion if 
one neglects arbitrary qnantities vanisking at least as c © 
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The dynamical variables describing the fluid are the Euler four-velocity field u(x") 
== S satisfying the condition ggg u^ uf = 1, the pressure p (4^), and the density 
o (x") satisfying the equation of continuity: 

(eur), = 0. 


We assume that the fluid is barotropic, o is therefore a given single-valued function 
of the pressure, and we assume that outside finite regions of space ọ = p = 0. The 


expression a (e i SAM p) describes the potential energy associated with the 
c 


compressibility of a unit volume of the fluid. 
The requirement that the action W be stationary under the variation 0g, leads 
to the field equations?) 


Re? _ 4 gP R= — 8x T”, (3.2) 


where 
d 1 
1 d 
re = Lg |o (ie å f uw or], (3.3) 
0 


If, on the other hand, we require that W be stationary with respect to a variation of 
the Euler variables in the class of motions satisfying the equation of continuity, we 
would obtain the equations of motion of the second kind for the fluid. 

We now come to the realization of the program formulated in the previous section. 
We introduce for this purpose new dynamical variables v* and u which are more con- 
venient in applying the approximation method 


ve 0 
= Ce a v^ DE , U = i , 
ce 


e= = 1 (ur va vh. (3.4) 


As may readily be seen, the density u satisfies the equation of continuity 


usa 


Ho + Qi, = 0. (3.5) 
The action (3.1) in the new variables becomes 
x 
u-&w= | a, (3.6) 


where 


Lm ly + es 


5) E the variation of w the motions satisfying the equation of continuity and the condition 
| Sap u% uP = 1 should be taken into account. The derivation of thetensor (3.3) from the action can be 
` found in the book of Fock $847 and 48 (Fock 1955) 
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L = — f dx u (gag u* uon 
1 3.7 
L= = | eyes, G0 
p 


EAE 


o 
In accordance with the approximation method we expand the field g,, into a series 
(2.9), while the dynamical variables are associated with the particular orders in c7; 
we assume that z is of the order 2, v? of the order 1, and p of the order zero. In order 
to write the action in the lowest order of the approximation method we must find hog 


and Ah, on the basis of the field equations. The corresponding calculations are given 


2 
in Appendix I. It is found that 
å i À , u s AE hor: 3.8 
hoo 2 f aw PSP 3 de Orn = ( J ) 
where | x —«' | = [(x* — 29) (x* — x4)]%, w = u (x), the integration being made 
over the entire three-dimensional space. 
After elementary calculations we obtain from (2.2) 


MN 1 
4 4 2 2 


In order to check condition (2.10) we should estimate to an accuracy of the fourth 
order the expression 
2V—e G 6* 
9 Gapim 
(in calculating this expression we make use of the fact that g,5,, = ^,5,, and 
hp = Ô* gap). Since hoojm behaves at infinity like O(r-?), and hon and "loo like 


848 rm Room O* hoo 
2 2 


O(r-}), then the above expression behaves like O(r-8); condition (2.10) is therefore 
satisfied in the present approximation. Making use of the field equations (1.5) we can 
transform expression (3.9) as follows: 


G = 45744 ho + Vantin > : 
4 2 2 


where 
Wm = — å Angles Ron 
s 9 "oom 2 9 
We still need for the calculation of L, the expansion of the radical 
-4 
1 1 
(gag v* u) = I — DX art hoo: 


i 
On the other hand, to give L, it is sufficient to know the relation between o and u in 


. + 4 
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2 er 
the lowest order, i.e. 0 ~ = ú and V—g= 1. As a result 


1 1 
pu faut Ef ance | annt. 
P 
x dpa 
0 
Jb s= I d. hoo + l dr n" doo 
fit ads 37.0550 167 ui 29. 


The surface integral ó W,,n"d,0 is taken over the surface at infinity. Since at infinity W, 
4 4 


behaves like O (r73) the contribution from this integral vanishes. Inserting into (3.10) 
hog from Eq. (3.8), we obtain the Fokker action for the fluid in the Newtonian order 
2 


0 
*(2) 


1 1 , be 
Up= f a . fax ecl [ar f ae PERS + 
0 
"D a (3.11) 
x dp at 
-áf wu |+ 2 Ja] 


We require that the above action be stationary with respect to the variation of the 
Euler variables 6x* vanishing at the ends of the interval of integration in the class of 
comparative motions satisfying the equation of continuity (3.5). Very natural 
inclusion of the equation of continuity as an auxiliary condition to the action principle 
can be obtained by passing from the Euler variables x* to the Lagrange variables 


a* in the description of the motion of the fluid: 


ar = (x9, a^), 
x? =a% for x9 = 0. 


In the Lagrange variables the equation of continuity has the well-known form 


(3.12) 


DE) „0 Ter 

Da " på (a, x (20 a*)) = (0, ab; (3.13) 
D(x). : S I 
Dia) is the Jacobian of the system of functions a* = x* (49, a’). 

We now change to the Lagrange variables in the action U p; simultaneously, we 
+3 
k 

restrict ourselves to motions satisfying (3.13) (ve also recall that v* = ca. 


Asa result we have 


gu [ao [- [oom f anne ee 5 
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(x9, x*(x9,a!) 


MEET fae fa daa - HS DALE) fa a a fl ce | 
|x (39, al) — x (x9, a’)| — ct pa ec) ; (x, xF(x?, a!) 


+3 da 5 LP (ata ay | 


If we require that the above action be invariant with respect to the variation 
6x*(2°, a*) we obtain 


© duk l 
fen. a’) 5 nn 


]«t 
e «fa ga u (0, a’) u (0, a”) = je Ook + 
0, a 
ae dat aah apy Pi (5, CA a) dt + (5.14) 
D (x x 


x) 

9x (x9, a 

5 | aa 6; 
0 
AD 


The third and fourth expressions can be simplified with the help of the equation of 


— 0. 


k 
continuity. The variation of the Jacobian ô p? can be calculated as follows: 
a 


D(x) KV på 8óx — D(x*)90x | 
D(d) "dar * Bak dx D(d)9w"' 


Ó 


S 


| x 
MF is algebraic complement of the element Dak of the Jacobian. In (3.14) we change 
over to the Euler variables and obtain 


xC) 


Jel- fene fe faw, uw = ghe Jp 
*() 
«2 
+ Å aspe] + fino =10, 
x) 
do 


= vo Tot v is the substantive derivative. The last two expressions above 


ý 
+ 


di? 
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vanish, and owing to the arbitrariness of the variation dx‘ we obtain: 


Bet 1 x 
— p (veg + ot v) + [ax un eae PG (3.15) 


These are the equations of motion of a perfect fluid. These equations together with 
the equations 
pp + (ur*), = 0 (3.15?) 
p = pl) 
constitute a complete set of equations determining the motion of the fluid in the Newto- 
nian approximation. 
We now come to the construction of the Fokker action for the fluid in the post- 
Newtonian approximation. After an elementary, but rather lengthy, calculation we 
obtain 


1 3 
. G = — — hom Room + hoo hoolm Roolm — > Polo Žoolo + 2 Rooim homlo + 
2 2 
> 6 2 2 2. 2 2 237172521 2 3 1 
1 


Í 
2 3 3 2 3 3 2 4 


If we apply (2.12) to calculate the integrand in condition (2.10) to an accuracy of the 
sixth order we obtain 
8y—g | 
VER Ó* gag = — hoolm O* hog + 2 hoo hoo|m O* hoo + 2 homlo Ó* hog + 
Saf|m 2 2 25403 2 3 2 
— Amir Ó hs + horm 6* hor — hoo|m Ó* hog — Room Ó* hog. (3.17) 
3 3 3 3 2 4 4 2 


Since for variations in Eq. (3.17) all coefficients behave at infinity at least like O(r-?) 
(see Appendix I) and the individual variations like O(r-1): then this entire expression 
behaves at infinity like O(r-?), that is, condition (2.10) is fulfilled also in the sixth 
order. 4 Ed: 

For further calculations it will be more convenient to use the expression which 


we obtain from (3.16) with the aid of the field equations (A. 5), (A. 14) and condi- 
tion (A. 18) 


G = 4r u ho + 82 u hoo + 87 u v" Nom — Ant u hoo hoo + 
4 3 247 


~6 2 
1 
+ g Pott homm + „mim + Walo: (3.18) 
where 
N + ! 
om "E oolm Aoo — oo|m "too — 3 hati hom + 2 hoo hoo olm =; 
1 
= hy (oti, — homi), 
Js hor ( tor om) 
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In order to find L, the radical (Eag v* v)? should be expanded to terms up to the 
fourth order ~§& 


Í 1 1 1 
z sæ pB\} = 1 — — 94 ya — — — — 
(gag v* uË) "V Y y ho +“ he 8 oo hoo 
3 
[8 (ve pa)? + a hog v? v* + hoa V2. (3.19) 
| 2 3 
Inserting into Eqs. (3.7) for L, and L, the expressions (3.18) and (3.19) we obtain 
-> 6 + à 
L Lo == d l d Ha 
ors een are: so FRAG US a Y dX B Boo + 


E san (0909 — x dix hu hog 3 fa 3% 4 hoo vt yt + 
1 
- f annie encom re lez fa % hort homm + 


1 
+ — Ir Le n" = er (far r) (3.20) 


The last term, being a time derivative, does not have any influence on the equations 
of motion; on the other hand, the integral $T, n" do = 0, since J, ~ O(r-?). In 
+68 - +6 


addition we see that in the above expression the terms containing ho drop out. Thus 


to calculate the post-Newtonian action it is sufficient to know Ayo Aa hom calculated 


in Appendix I: 2 


ee ee ss , wm 
hoo = far gs hon if: Eum RN S. 


and to know only the asymptotic behaviour of ho) necessary to estimate the expression 


4 
(3.17) and W,,. Because of this the method of deriving the equations of motion by 


constructing the Fokker action is considerably simpler than previously used methods. 
In 1938 to find the post-Newtonian equations it was necessary to know oo hon» oo Pm 


and hom SR at present the first two of the above expressions (of Nm ESI 


TUM than the others) allow one to determine the equations of motion. 


We determine L, also from Eqs. (3.7). For this purpose we must know the rela- 
+6 
tion (3.4) between o and u with an accuracy to the fourth order: 


— e 1 1 
pe dices (i- Lenis) 
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and 
| = il hov: 


Next, we determine the Lagrangian OG L, + L, + Lp insert into it the values (3.21) 
EE 
of the fields hj, and Aom; we EG the integral f hoy; homm dax by using the 
2 3 Sind 


identity 
da x TA IRR = —27 |x’ — x" | + constant 
|x — x| |z — x" LD 


the fixed constant (= oc) appears in L in the quantity to be differentiated, and therefore 
+6 


does not make any contribution to L. Ås a result we obtain 
-6 


1 
DE dzu +> d4x u VI v? + JE 


3 uw ; 
d ana la pa a ra 

TRA NOH ze + v'« v'a) 2 fas [ar m vem xe + 

1 d. x! OT REN ” L4 1 , uw 

F | dar Jar mo” + > | dax | dgx Re: + 

1 1 
——[4x][a Aue, u 

5 får far far wee Apa TEST) | kle" 


1 1 z 
x Lå 20 — x p 
TRE ini ne : Ge ei 


c? 1 , 
LISE Tide pal , u 
T NER Jesi 


The interpretation of the individual terms in the Lagrangian is rather simple. The 
first expression represents the rest energy, the second and third represent the kinetic 
energy together with the first relativistic correction, the next three terms describe 
the post-Newtonian gravitational retarding effects, the seventh term is the Newtonian 
potential energy, the eighth term is also of the potential energy type; to obtain here 


the interaction between points x and x’ the elementary interaction between three 
pl 4} 


E u ; à š 
kr should be integrated with respect to x” over all space 


where 


points 


k m 
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(here there is a certain analogy to the interpretation of quantum transitions 
in the quantum mechanical perturbation theory as a complex transition through virtual 
states); this term was given in symmetrical form with respect to the variables x, x’, 
and x". The last two expressions in the Lagrangian (3.22) describe the previously 
mentioned energy associated with the compressibility of the fluid. 


If we require that the functional Up = f dx? L be stationary with respect to the 
+7 +6 
variation 6x* for the class of functions satisfying the equation of continuity (3.5), 


we obtain the post-Newtonian equations of a perfect fluid 


P 

Le ,_ M 1 (de), DEUM. 

[e (14 pore ts fae en D) aif ae Eere FF 
0 

er dax' |x — x' luy un’ v's +l lg-L ww 43 po d jJ 

RT EL us "RW * 2 ° Ix—x| 

P 
1 dp uw " 1 , , ler 
en 2) oa f age HH ont da" |x — "lige mye! v5 | rti + 
0 


1 , 3 , t 
+ E 3% uw 1x — x |ipay v? v — | dax' = T ap 


3 ies 
bee tat La fie n (ren et 


u 
pe eee re wre 2 [ar za 


x—x|ix — 2 
EG | AE d) +2 (w + up} | =0 (3.23) 
EN p es 
0 0 Š i 


(the calculations leading to these equations are similar to the calculations in the Newto- 


nian case). 
The full set of equations determining the post-Newtonian motion of the fluid 
is obtained by supplementing equations (3.23) by the equations 


Po ot (uv?) = Ü 

p = pip) (3.23) 

We note that equations (3.23) have a different content than equations (3.15) 
which appear to be the same. The dynamical quantities u and 2° appearing in them 
are in both cases defined to an accuracy of different orders. In Eqs. (3.15) v* satisfies 
the Newtonian equations of motion, and in (3.23) the post-Newtonian equations of 
motion. On the other hand, the pressure p occurring in the post-Newtonian equations 
of motion is only the Newtonian pressure. This can be shown most easily with the aid 
of the Lagrangian GRE Let us assume that the pressure p — 3 er + ..., where p 
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is the Newtonian pressure and p is the post-Newtonian correction. p can enter into 
2 2 


the Lagrangian only through the expressions — > JE 3% U få — and — a fa gt D. 


Let us consider the first of these expressions; since 


T s pt 5 


fz- j m : fes od dp _ festa, p x å 409), 
o PG 


then 


Hence the post-Newtonian correction to the above expression cancels out with the 
x Li * ` 
correction from T ifi d,x p, that is, the pressure appearing in the Lagrangian (3.22) 


is the Newtonian pressure. 

The set of equations (3.23), (3.23’) could find use in the case of fast motions of 
fluids in a gravitational field; it would seem that that these equations could be used 
e.g. for fast motions of fluid stellar material, or in a hydrodynamical description of 
the evolution of super-nova stars shortly after their explosion, when the gravitational 
forces are responsible for the mechanism of their evolution and their velocities are 
still large. 


4. Drop model of a system of n bodies. 


The equations of motion of n point bodies described by singularities in the ten- 
sor T*? contain infinite terms (Infeld 1954) which can be removed by renormalization 
of the mass or by introducing into the energy-momentum tensor the so-called “good“ 
6 — functions (Infeld and Plebanski, 1956, 1957; Infeld 1957) whose mathematical 


formalism automatically removes the infinite terms). 


8) We notet hat the taking of „good“ 6 is not sufficient to remove the infinities. Here we need an 
additional formal procedure. In the equations of motion there appear expressions of the type 


2% T T 
|x — x'| WESTIN 


whose limit for £ — 0 is infinite for each model of the ó function. In order to remove the infinite terms 
one should introduce the independent parametrized models of the å function Š, (4) and Š,, (x’), then 
one can choose the models ó, (z) and 6,, (z) such that 


AL ef ae OO mg 
ae 80 rer == Ed 
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The formal renormalization and the procedure of *good* ó however, require justifica- 
tion. It is also necessary to explain the physical source of the infinite terms in the 
equations of motion of bodies described by singularities of the field. 

To explain the physical side of the process of renormalizing the equations of 
motion, it is convenient to employ a drop model of a system of n bodies. In this model 
we assume that a perfect fluid described by the equations given in the previous section 
is separated into a group of n drops whose diameters are quantities of an order of 
a characteristic length /, where the length I is small in comparison with the mean 
distance between bodies. The mass of one drop denoted by the index A (A = 1, 2, ... n) 
is equal to 


A 
u = J dyx n (4.1) 
Q 


A 
where Q is the space region occupied by the fluid forming the Ath drop. From the 
A 


d A 
equation of continuity (3.5) it follows that X RN 0, and therefore the mass u is 


dx? 
A 
constant. The world line &^(x9) of the centre of mass of the Ath drop is defined as 


AA A 
UE) = f dysux*; EY) = 20. (4.3) 
J | 


Q 


Next we assume that the velocity field of the fluid in the Ath drop is homogenous and 
A 


equal to the velocity of the centre of mass £^. This assumption means that we neglect 
rotations and deformations of the drops. The case of rotating bodies is discussed in 


à separate paper (R. Michalska, in preparation). We assume here, indeed, that the 
A 


density and pressure of the fluid vanish outside the regions Q. 

Making use of the above assumptions, we derive the conditions of equilibrium 
of the drops for our further considerations. For this purpose we write the Newtonian 
equations of motion (3.15) of the fluid forming the Ath drop (we take into account 

A 


the fact that x € Q) 


new Y fa ut weer (at 


BFA B 
a 


+ Ir aX u(x) u(x’) = =) = | Pie = 0. 


The expression eu on the left-hand side of the above equations depend on the 
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A 
quantity I characterizing the dimensions of the drop. Since the mass u = m dx u 
À 
Q 
does not depend on I, then the density u, as the dimensions of the drop tend to zero 
(1 — 0), behaves like £3. If, therefore, we expand the left-hand side of equations (4.3) 
into a power series with respect to l, then the first two expressions in (4.3) will have 
an expansion beginning with 173. On the other hand, the third term in (4.3) is a quan- 


A 
tity of the order /-5, since for x and x’ in the region 2 the expression (s 1.) ip is of 
A B 
the order I? (for x € Q and x' € Q, BA A, this expression would be of the order 19). 
We thus conclude that pj, must be a quantity of the order 179, since otherwise we 
could not satisfy Eqs. (4.3) in the order 15. The left-hand side of Eqs. (4.3) therefore 
has an expansion in I in the form a. 5/75 + a. 4173 + ..., where each of the coefficients 
of this expansion, by Eqs. (4.3), must vanish. From a_, = 0 we obtain 


fax (x) u(x’) pec = E Pik- (4.4) 
p ipen 


a 


The above equations are known from the hydrostatics of the equations of equilib- 
rium for a fluid drop in the field of its own gravitational forces. These equations are, 
as it follows from their derivation, on the one hand a relation of the order c4 sa- 
tisfied to an accuracy of O(c~*) (O(x) is the symbol of a quantity which tends to zero 
as rapidly as x), and, on the other hand, a relation of the order 15 satisfied to an accu- 
racy of O(13). 

From the equations of equilibrium there result two equations which are important 
for further applications 


3x 1 ‚ulx) u (x ^ 
= dp = fax [ar — (4.4) 
4 4 A vod 
Q 2 a 
and 
ES P 
ps dp == 4 ux) Hn 
c u(x) = fa ve š (4.4 ) 
Teo dem 


The derivation of these relations is given in Appendix II. Equations (4.4’) and (4.4), 
as may be seen from the derivations, are relations of the order I satisfied to an accu- 
racy of O(i). 

Let us now return to the post-Newtonian motion of drops of a perfect fluid sat- 
isfying the assumptions formulated at the beginning of this section. Since the fluid 


A | 
occupies disjoint regions of space Q (4 = 1,2, ... n), the the volume integrals in the 
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A 
Lagrangian (3.22) are equal to the sums of the integrals over the regions Q, or 


1 1 
RU TRAE PAT TÉ 
ASA A A a 4 
Q 
3 
n [ef TU) EG PRE ZX Je ap EGG Be vay'a 
BA 
Q 


4 


1 
yy fel 3% uw x — x |y veut + T ADD IE x at 
Q 2 


= 


APE 
Q 


SD |e few [e (em i 
ed dre EE C 
6 E 3 3 3* J 3* lx— x" | lx — x"| ES ex] jx’ —2"| ai 
Q Q 


2 


P 
1 f å d, X 
+a yy a far 2 [yt fan 
4 B 5 A Ae dp 4 À 
Q Q 
2x 7? Hp 
B. a yy [a dax we: _ (4.5) 


(we retain here the notation u’ = u(x’), u” = u(x")). We expand this Lagrangian into 
a power series in / and write the explicit expressions up to the order /? inclusively. 
We transform the first three terms in (4.5) by making use of (4.1), and (4.2), and the 
assumption of the absence of internal motion in the drops. For BE. next terms we also 
employ the following expansions: 


(a) LE [se far t a So pA D i: Ju dax’ pes Ë 
Q Q 
(b) DE f dax | de up. |x — x lag = Nr uu 7 js +00) + (46) 
A B 4 B 


A BFA 


+ yy. Å dyx il dx up" |x — x jav 
FL tå 4 A 
a Q 
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« , tt , tt 1 1 
SEXE fus [ar [2: piede re a 
ARTE SCA B C 


Q 
ABC 
t kex] E =)= Då på på Hit (me + ape T din SF 
À B+AC+B ji MT fees 
C+A (4.6) 


EE Hs bats DE a (Je fargi 


A B+A r 
A u” 
th fax fax LH s Efefor für meer uti 
B B EG 
Q Q 
where 
ABI ACH AB bes 
> 
r= |; ¿l rje= ar 
2 £a 


In addition we introduce the notation 
A A 7 
u? o = fax [ax = =s (4.7) 
Á A & % 
Q Q 
A 


A 4 
(The expression u? w is a quantity of the order /-1, since for x and x’ belonging to (2 we 
have |x — x'| ~l) and we make use of the relation 


24,4 
A dx & — lier up = — 3 på O dap. 


Ås a result of the transformation described above L takes the form 


+6 
2rd ess 1 
L = Lu) + > D Ep £% (2642 [den f2)+ 
A 0 
Q 
1 A B B ^ 
1 1 : d 2 i 
225» (moth faxu [24% faxs)+ 
A B+A r ; B ° 8 
Q Q 
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dd 
Q o Q 
l å pM 3m dp 
UE 3% U =D dax | dyx Er pons 
A A A A A 
Q Q Q 
x 2x 'p 
+55 <s žy far [as z Z t 00, (4.8) 
4 A 4 A » 
Q Q Q 


where 


1 A B AB A 
—i» » Hu T\ £0 À Elo Bh + 


A BFA 
1 Å 1 ABC 
uu l | 1 
EO A 5 D nb (mt mm + ace) — 
A B+A r A Bt A CHB Der ndr ror 
1 AB AG BL) TUS 
— 3 Y weet) (as) (4.9) 
A B+C r 


is an expression independent of /, the expressions in (4.8) are explicit expressions 
of the orders I? and 17? and therefore are singular for / — 0, while O(I) tend to zero 


for 1 —0. 


We now want to show that the Lagrangian (4.8) can also be transformed in 
such a way that the singular terms enter into it only through an additative 
A A A 
correction ôu to the mass u in the Lagrangian Tees (a), where the lowest term 
A 


in the expansion du according to the EIH method is of the post-Newtonian order, 
MU A A 
i. e. ĝu = en oh gu +... We thus want to show that the singular terms in L can 
-8 
he AERA, d à renormalization of the mass. 
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A 
In order to determine the correction öu we consider the Newtonian Lagrangian 


A 
Leg (u) 
24 


A B 
A A 1 AA A il uu 
Lelu) = — Yi n 3 2, B ën "odo Y "AB ° 
> A A 
A 
and in it we replace u by 
A A A 
reme Ôu, (4.10) 
4 


A 
where du is a quantity of the order ct. We then obtain 
4 


A A A A 
Lg (u + du) = Lag (u) = 31 58 + O (7) (4.11) 
24 4 +4 A 4 


reg 
4 


we obtain by discarding in (4.9) the post-Newtonian terms) 
A 1 A A *d 
pom ES 2 Ed 2.5 ce 
Len + Oa yf dr å [axe f. 19 
^" T å ö 
Q Q 
when we take 


b 
A ] 4-4 x X dp 
4 0 


The renormalization procedure given here for the Newtonian Lagrangian would not 
have any meaning if we wanted to stop at the Newtonian approximation. For, terms 
of the order I appearing in (4.12) do not depend on time and have no influence on 
the motion of the drops, and therefore on the motion of the material points into which 
our drops pass as I > 0. The correction (4.13), which we determined when considering 
the Newtonian Lagrangian, however, now permits a meaningful renormalization of 
the post-Newtonian Lagrangian. For, we have 


A d rA 
We see that Leg (u + ôu) will be identical to the Newtonian Lagrangian L (which 
> 4 +4 


A PET. A DET OS nt Met 
FAD = Treg (4 + óp) + 2 deg u Ei €fo + 


1 1 b B rsi BO JA 
N N as won + udn), | (4.14) 


A BA r 4 


The General Fokker Action Principle 355 
Inserting (4.14) into (4.8) and making use of (4.13), we obtain 
A A 1 AMA ] 414087 
= ls Ó = a Et) [— — u? — 
= Des (u + 1 +5 D | go à fax) + 
Q 
1 1 [4 LIB NB 134 | < 
HA LP o [aeo +e (75 pu DLE 
4 BA T ` 
Q 
3x 
ve j ten |- Y 64 +09 +00, (4.15) 
p A 
Q 
where 
hes l 23% A eoi : 
Pow ed EHE à 
À À 
Q Q 2 
r : 
1 v qm E 2x | ap 
a fees fae katina SK 
À B 0 4 A 
Q Q Q Q 


is a constant of ihe order 17. If in (4.15) we make use of Eq. (4.4') satisfied to an accu- 
racy of O(l) then we obtain, by (4.7), 


A A A 
+6 +6 Are 
4 
The quantity 2 lp as a constant can be omitted from the Lagrangian. If we now let 
A A A 


the drop eom go to zero and assume that the physical mass uso, = U d % 
of a point in the gravitational field is finite, we obtain 


A 
lim (L 4- 2 on) = Lees (Bord: 


1-0 +6 


ie. the well-known post-Newtonian Lagrangian of a system of m material points 


represented by the singularities of the field in general relativity theory. 
A B 
We already determined the constant dm which normalizes the Lagrangian in the 
4 
post-Newtonian order on the Newtonian level. By analogy with the mechanism of 
A A 
obtaining du we should expect that the constant dy will be a renormalization constant 
4 6 
for the post-post-Newtonian Lagrangian L. We can be consistent with the procedure 
28 A 


presented here and not leave the constant 2 ou out of the Lagrangian L, but include 
+6 
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YR A 
it into the Lagrangian Leg (u + ôu) as a further correction to the mass. By (4.9) 
+6 4 


we can write Eq. (4.17) also in the form 


A A A ; 
L = Ly, (4 + Ou + 6p) + OG) (4.17) 
+6 —+6 4 6 


(since where the difference between the right-hand sides of Eq. (4.17) and (4.17) 


is of the order O(c9)), and we assume that 


A A Å A 
Unorm. = H + i. A d 


Since, as we have already noted in Section 3, the pressure occurring in the Lagran- 


A A 
gian L, and therefore in the constants ôu and dy, is a Newtonian pressure, then, by 
-6 4 6 


employing Eq. (4.4) — (4.4) we can represent the renormalization constant in the 
form of expressions of the potential self-energy type 


43:2] r he 
= + [ax fax PET (4.18) 
A A 
Q Q 
4 1 , H hows i Eg | 
|, Jar dax" uuu nat: (= £l: Å 
Lo] 


We note that if we had wished to describe the motion of a system of material 
points in general relativity theory as the limiting case of the motion of drops of a conti- 
nuous medium whose pressure vanishes everywhere (p — 0), then: firstly, we could 
not satisfy the equations of equilibrium (4.4) and therefore the Newtoniane quations 
in the order l-5 for such drops, and secondly, we could not renormalize the post- 
Newtonian Lagrangian, since it would not be possible to bring together the terms 

A 


of the order 17" in the expressions entering the Lagrangian Lies (u) only as corrections 


A 
bu to the mass u. 


* In conclusion, the authors express their gratitude to Professor L. Infeld for the 
kind interest he has in this work during its preparation. 
Appendix I 


We shall use the EIH approximation method to solve the field equations (3.2) 
which we now write in the form 


CRE = —8z (re > - eh r) , (A.1) 
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where 


R£ — V—g Red — V—g g^" gå he (A.2) 


Me dvd 


Inserting into the above g.g = Nag + hag and expanding h,, into the series (2.9) 


we obtain the expansion R% (making use of g,, g” = 65, we express in the above 
expansion A% by hap)» Where in the lowest order of the expansion we obtain 


Il 

R»—-———h rr} 

9 z 0! 

Rm= lp Sende mom h hd: (A.3 

Å 2 9 Oolmn 2 Men Fi) nier As > . ) 
1 

Ron — — 9 (hom|rr ae. hor|mr SE hrr|mo vas firmior). 

3 3 3 2 2 


In expanding into a series the expression 870* = 8x (T'*? — 1 g*? T) on the 
right-hand side of Eqs. (A.l) we regard the dynamical quantities as quantities of 
a fixed order in c}, and we expand only the field h,,. If we assume that p is a quantity 
of the order zero, v? of the order 1 and wu of the order 2, we see that the lowest terms 


in the expansion of the tensor T* in (3.3) are T'?, T?" and T””. Thus the expansion 
202 4 
of the individual components 0 begins with 


1 1 

LT oa ae 
6) p z Mi 
Com qom IE (A.4) 
3 3 

1 1 

mni Sn 00 —— a n 

e 5 Ó n 5 om" y, 


Having now the lowest terms of the expansion (A.3) and (A.4) we can proceed 
to discuss the field equations in the lowest order of the approximation. For «f = 00 
we obtain 
hg, = 87 T9, (A.5) 
2 


ie. Poisson’s equation. Its solution is 


d 
E E 9) d ! poo ; 
hoo fas [x — x'| 
; (A.6) 
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where 
qos = fu Es 
2 2 
M = uix). 


To particular solution (A.6) we can add the general solution of the homogeneous 
equation hoss = 0, and therefore an arbitrary harmonic function. The solution of 
2 


equation (A.5) must, however, satisfy condition (2.6) (since y vanishes outside finite 

regions of space, then solution (A.6) satisfies this condition), and therefore vanishes 

at infinity, while only the non-singular harmonic function vanishing at infinity is 

zero; on the other hand, we discard the singular harmonic function for physical reasons. 
For af = mn in the second order we have 


ll 
xx d (Ayrimn == limner Renta EE been 3 = hoolmn > AO mn E» (A.7) 
2 2 2 2 2 2 2 


By means of (A.5) we can readily check that the particular solution of the equa- 
tion (A.7) is 


Ran = Dn hoo 
2 2 


mod Hi egre Eden (A.8) 
2 


lx — x'| 


Pre = 
* Dm f as |x — "| ° 


To the above solution one can add the general solution of the homogeneous equa- 
tion corresponding to (A.7). This homogenous equation has the form 


(nr = Anim) |r FG (^u, "a Pip) n == (A.9) 
2 2 2 2 
Equivalent to the above is the equation 
Aan)r xn him — Jmr|n ° (A.10) 
where fm, are arbitrary functions of x*, where from it follows that fn, = —f,,- From 


(A.10) we obtain 
Late i rete Sa a = 0, 


or 
(E 73 br d byr 


(bm are three arbitrary functions of x°). Inserting the above equations into (A.10) 
we obtain 


Pennie En bn =F LE ER 
2 
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A Flur E IS ea 


or, using the notation 


we have 


S, S, 


mn|r TT rmm T 0. 
From the above we obtain $, = q,j and from the condition of symmetry Spn 
Em Guns nat follows that 'g, — Pim (p is an arbitrary function of x). 


Hece we can write (A.8) in the form 
h. = Dan x Dei e Pinm = amin fe Om , (A.12) 


where a,, = bm + $ |, are three arbitrary functions of x*. From the above derivation 
it follows that (A.12) is equivalent to the equation (A.9), and therefore represents its 
general solution. The general solution of equation (A.7) is therefore the sum of solu- 
tions (A.8) and (A.12). From the third of Eqs. (2.8) it follows, however, that there 


always exists a transformation of the coordinate system x” = x” — a"(x*) which 
2 


will reduce the general solution of Eq. (A.7) to the form (A.8). In order to restrict 

the class of coordinate systems introduced in Section 2 we impose a condition on the 

generator a”: Amn + anm = 0; we thus restrict ourselves to the class of reference 
a # 2 


frames in which the general solution of Eq. (A.7) has the form I = 0, og: We 
2 2 


still have at our disporal, of course, the free of the generators ag, dy, Ag, ..., Aps Que 
3 D G 4 5 


We determine h,,, from the equations 


I 
9 (Aomrr SE horlmr m hrrimo = Arm|or) = —8x T°. (A.13) 
3 3 2 2 3 
Since h,, = On Moo, we then obtain 
2 2 x 
LA TE Por ze 2hogjom = — 16T". (A.14) 
3 3 2 3 


As in the case of Eqs. (A.5) and (A.7) the above equation is already not integrable 
for every motion. Å necessary condition for the integrability of Eq. (A.14) is the 
following relation between the dynamical variables: 


hoo'omm + 82 Tn = 0. (A.15) 
2 3 


which is obtained by differentiating both sides of (A.14) with respect to x”. By (A.5), 
Eq. (A.15) can be written 
T^ + Tw 10. (A.16) 


The condition of integrability of Eqs. (A.14) is the satisfying of the equation of conti- 
nuity in the Newtonian approximation, 
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On may readily verify that a particular solution of Eq. (A.14) is 


1 uom 
= qom = '—P A.17 
fs s fart lua] afas xx ae 
By (A.16), this solution satisfies the condition 
Zhoojo — homm = 0. (A.18) 
The general solution of the homogeneous equation omirr — hos, = 0 corresponding 
3 


to Eq. (A.14) is Colm where do is an arbitrary ede. of the variables 4^. Thus the 


general solution of Eqs. (A. 14) has the form 


, 


EUM ESN A.19 
30 f Ce mx] | en ( ) 

Comparing the above with (2.8) we see that (A.19) represents the solution (A.17) 
of Eqs. (A.14) in a coordinate system which can be obtained by an arbitrary transforma- 


tion of the form x" = x? + a9. This arbitrariness is not required to narrow down 
3 
the allowable reference frames by further restrictions, since the post-Newtonian 


equations of motion are invariants under the transformation x'? = x? + a?. In order 


3 
to show this let us investigate how the post-Newtonian Lagrangian L = L, + L, + L, 
> +6 +6 +6 
transforms. By (2.8) the above transformation transforms the field h,, as follows: 


The Lagrangian L,, as may be seen from (3.7) and (3.16) depends on the fields hom 
26 3 
and A, through the intermediary of the combination 
4 


Room homo pz Room) a ery at ==, homi) (Pom TE horm)» 


or Le — L, L, is a function of ho only, and therefore it is also not subject to the 
= id nd. 


transformation. On the other end L, (see 3.7 and 3.19) is a function of hom and hoo 
through 4 hoo + Road“ , and ee på = Ly — ojo Hence L’ = L— ojo; as 
E 6 6 -6 
a result of ui DUREE generated by a, the ul time derivative ojo is added 
3 


to the post-Newtonian Lagrangian, which has no effect on the aptas of motion. 
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We shall now discuss the behaviour at infinity of the quantities hyp and their 


derivatives. In our case T° = y and T?" = uv”, where u = 0 outside the finite region 
3 


2 
occupied by the fluid. From Eqs. (A.6) and (A.17) it is seen that for | x | = (x? x?) 
considerably greater than the dimensions of the region occupied by the fluid 


x— x [3 = (x)! + O( | x |-*), and hence hog and À, behave at infinity like | x |-!, 
= et = Ae SUM 22 


and Ag, and Roms like | x |-? (we shall indicate this briefly by hoy ~ hom ~ @( p 
2 3 2 3 


Om|0 
2 


we must make use of the Newtonian equations of motion 


Room ~ Roms: ~ O( x [-?). In order to estimate the asymptotic behaviour of À 
3 3 


T^, — Tn, + T bo + + hom T” = 0. (A.20) 
+4 4 3 2 2 


For, we have 


1 1 
= P EEG E "IT |km 
[art 3% Lx — «| la Å == > 
+4 I dx | dx” rad al (eas icd und 90 T0 
ç Rx) N xz oz 


(here we took into account the fact that T^" = 0 on the limiting surface and we made 


4 
use of (A.6)). For | x | large in comparison with | x | 


1 1 ji ; på x i 
m — "Ten. +O ed LE of Ji 
homo Tf omm får fax [r x ) i uo 


Thus it follows from (A.20) that homo ~ O (|x|?), since 
3 


J d, a J d, x" m T'® T1100 — 0. (A.21) 
24 9 š 


[x VEG x 
We did not need to know Ag to derive the equations of motion. For a discussion 


of the existence of the Fokker action for the post-Newtonian equations, however, 
we must know how the quantitites hog, hop and Ó*Ag, behave for large |x|. The 
4 4 4 


behaviour of the expression WV, in Eqs. (3.18) and (3.20) also depends on the asymptotic 
ET: 
behaviour of Ag. 


4 
ho satisfies the equation 


4 
/ R” — — 87 9%, 
4 4 
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where x 
R” = V—g g™ g% Ry 
š 4 

By the last of Eqs. (A.2) we obtain 


il 1 
Ry = 2 hols = > hmnioo + homlom T Aco, 


where A, is a nonlinear expression constructed from one of the expressions of the 


4 
type hnip Pars and Ayn Apgirs contracted in different ways for all indices. Hence 
2 2 2002 

1 


1 , 
R» = — 2 hoo|ss == ° hmnloo zs hm|om + A OG 
4 4 2 3 4 


where A’ is a nonlinear expression of the same type as Aoo Ayo thus satisfies the 
4 4 
following Poisson equation: 


Roots misa Run|oo “is hom|om +240 +167 0%, 
4 2 3 4 4 


or, by (A.18), 
1 ; 
Pools Fr D homom AT 2 A 2 = 16 n (A.22) 


© vanishes outside the finite regions of space. The contribution from @% to hoo 
4 4 


4 
therefore behaves at infinity like | x |!. Expressions of the type À, Magis in Å also 
2 4 


palrs 
2 


give a contribution of the order | x |-!, since after the contractions and after taking 


into account (A.8) they are transformed into the terms hog Agoyn = 87 hoo T”, which 
2 2 212 
vanish outside finite regions of space. On the other hand, expresssions of the type 


h mnp hairs in A’® behave like Lg and their contribution will be of the order 
2 2 4 


|x|-?. We have greater difficulty in estimating the contribution Ao to hog, in (A.22), 
4 4 
coming from ho — | x |"? since there is the danger that this contribution will 


behave like |x| In |x|. We shall show that such is not the case. The dynamical 
variables entering into Aj, satisfy the Newtonian equations of motion (A.20), which: 
3 


leads to the result that also hoo ~O (|x[-?). Since ent = $ Romom” then 
4 4 3 


= 1 , [A 
ho=— far V omi 


from which, by (A.17), we have 


= 2 1 1 
TEEN, d , d Ht T" 9m Oe ster kø ie ENT 
hon zf = J "eie" Fa. lx— x'| . 


x — «|. 
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With the help of the identity 


I 1 
fax te =; Ex] = 27 E: — &" mr 
we obtain 


he = — 4 [ da Tes e—a" me — à f dg Tn a | py + 
un Hi dx" Td |x — x" line SA D f dax" iid 54 ka ^00 ix — x" luy 7 
In transforming the second of the above integrals we note that T”*” = 0 at the 


4 
boundary of the region of integration; in the third integral we use (A.6); as a result 
we have 


Tee EN Hn pu " ” T ” 
lema [an Dena lå | pr | dye Tile — 
+4 4 


tt ur , il 
1 fay lx — x lec faer (m, 3 


For |x| large in comparison with |x'|: |x — x'| ~ |x 


xh xk ; 
" cl) Ix — x | Im’ 
TA + O (k9; le — hum ~ O (2-3), therefore 


EB 
im ferme Jae fa serm rør et ol) 
from which, by TS and TR it is seen that hog — Q ( (Ix|-3), and hence also 
oo ~ Q (|x|3). On the other hand Ago, behaves like Q (xf?). The variation Öfh,, 
behaves at infinity, of course, at ER like the individual fields h,,. Hence it follows 
from the above that in formula (3.17) a y—&G ö*g,, ~ O (|x|) and in (3.20) 
xus r. SN ap m 


Appendix II 
We now derive Eqs. (4.4) and (4.4") from the equation of equilibrium (4.4) 


J dx p (x) u (z) = |& == = Pik (B.1) 
4 2% 
Q 


A A 
We multiply both sides of Eq. (B.1) by (x — &*) and integrate over the region (2. 
Owing to the obvious identity 


j A A 
(at — & py = [G5 — ë) ply — 3p. 
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4 
and owing to the fact that p = 0 at the boundary of Q we have 


pm JEG u(x) u(x’) (a a)r (å — By 


(B.2) 
k ^k y 
xi Eid uc) w(t) EIES = > 
since x* — B (24 — x) + P5: — Ë, then the integral 
xk — o/h aS: 
pal ART ue) ue) EEG ë | 
63 


EE 


pho ax a a får får ux) u(x) (= É LUE CEOE 


In the last integral on the Sai land ks we interchange the roles of the variables 


x and x'; as a result 
e heel: ses LIU + 
- AA rp ibi 


or i E 2 1 Li hac HAV SA 2 nr a "Y 


Fon på 03. Ai (84) ve oe t men sad a odi. id s 
` Y ue ip å - ae i £ à y å "u : 
| _ +] eria | 


ALI) 
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FIRST FORBIDDEN BETA TRANSITIONS FOR THE UNIVERSAL 
FERMI INTERACTION 


By Bronıszaw KucHOWICZ 
( Received February 16, 1959) 


Formulae are presented for the electron polarization in recoil experiments for the 
case of first forbidden beta transitions. The formulae have been derived by means of the 
Cutkosky-Deloff trace technique. All Coulomb effects (for a point nucleus) are included. 
In accordance with the present knowledge of B-decay the vector and axial vector parts 
of the universal Fermi interaction are retained, with different coupling constants. 

The formulae for the case of 0 — 0 (yes) transitions are compared with the formulae 
given by Berestetsky et al. Some discrepancies are shown to exist in the asymptotic form. 
An attempt has been made to explain this. 


Increasing experimental data now tend to support the validity of the universal 
Fermi interaction which was proposed by Feynman and Gell-Mann (1958) and Su- 
darshan and Marshak (1958). Many experiments have been reexamined since that 
time and the theory has found a strong support in the B-decay and in the leptonic 
decay modes of the mesons. The theory has also provided new interesting effects 
(e.g. the *weak magnetism“ of Gell-Mann (1958)). The possibilities of further applica- 
tions have not been exhausted, even though some of the effects may be very weak. 

The electron-neutrino angular correlation connected with the measurement of 
electron polarization for first forbidden beta transitions seems to be a problem which 
was hitherto only partially solved. In addition to references 4—8 cited in a previous 
paper (Kuchowicz 1959) one may mention the paper of Berestetsky et al. (19583) 
dealing with unique transitions, Coulomb transitions with AJ£.2 and 0 — 0 (yes) 
transitions, as well as a previous paper of the author (Kuchowicz 1959), where unique 
transitions are treated. The problem is fully solved only for unique and 0 — 0 (yes) 
transitions. One also needs, however, formulae for other first forbidden transitions 
with AJ = + 1,0. 

In the present paper general expressions for electron polarization in first forbidden 
beta transitions are derived. In accordance with the present knowledge of B-interaction 
only the vector and axial vector couplings are retained. The universal Fermi interaction 
is, however, generalized in such a way that the mixing ratio of the V and A couplings 
may differ from unity (and may even be complex). Therefore the two complex coupling 
constants Cy and C, are used throughout this paper, while it is assumed that the 


(547) 


348 B. Kuchowicz 


parity non-conserving coupling constants Cyay and C42, (in the notation used in 
the author's previous paper) are equal to Cy, and C!) 

The trace method of calculation is employed and the correct Coulomb functions are 
presented as the product of a certain operator Å times the usual plane wave spinor. 
The notation and the coordinate system are the same as in the author's previous 


paper. 
The squared transition matrix element from an initial nuclear state to a final 
nuclear state with the emission of an electron and an antineutrino has the form: 


| < H(p, 4, Z, Ü > |? = F(Z, E) [V + nN + t, T, sind + t, T,sinó] () 


The longitudinal polarization of the electron is given by 


P, = 7 (2) 
while two transyerse polarization components are: 
d Las: : 
PE W srg and per 7 Sin 90. (2°) 
The functions W, N, T, » and T, are given below: 
ae 4 Ma 
= [Ga 247; Aut [Cr]? D yM; Vag + 2Re (C Cy) 2; ave 1; + 
- i=1 i=1 
+ 21m (C4Cy) Sat, T ati 


= Gal: Lam mai ue Vag + 2Re (C45) Bay By Ta + 


+ 2Im (C267) > Av MN, d To ati 


1) A generalization of the formulae for the case in which zy, and «4 differ from unity is possible 
if one makes the substitutions 


1C;2?#>4${¢, PQ + | ; |9) 
C; Gj > 1 C; G (1 + a; a) 
in the expressions for W and Tz, and 
| C; |? 1 | C; |? (z; + a) 
1 C; C; > 1 C; G (z; + a) 


in the expressions for N and T,. 
The subscripts denote the kind of interaction (4 or V). 
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6 4 


4 
RF cyu sts RANG NR 
i=1 ici ; 


i=1 
å 
+ 21m (CaCy) > av; a: T3443 
1=1 . 


6 4 


4 
D = CA > aA, 1% > RM Ag, 2Re GK > el 
aem 


i=1 i=1 
m t 
+ 2Im (C4Cy) p Av Å; Lua 
i=1 
The M; are expressed by means of the nuclear matrix elements in standard notation: 
AM = | f Ys |? 


es 


48 =| f ox?| 

Jn = >| f Bal 

a, =| fap 

vr = | f Ff 

vM = 2Re [(f e) (f *)] 
em = 2m (JJ) 


DEED) 
No special assumptions about the relative phases of nuclear matrix elements have 


been made. If time- reversal invariance is asumed to be valid for nuclear forces, the 
interference terms 49ft 9t, , ap, and sy Mare to be put equal to zero. The coeffi- 
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cients Ap; Vp; and J,; are listed below: 


Au = 2 (Lo + Ly cos 9) 


4 
As SIS —39M + 2M, +3 2 — 4 a Ni + 218] cos Ÿ 


2 
Toye 
As = —2Ro cos 3 
2 2 , , 

A= z qL e —2N + |2 7-22] cos Ÿ 
pote free tee i : Š 

1699 d ot q Ka,-2 3 — cos å + s Mo + 2My + L + 

1 QUEE. Bem 2 I 

*[- 3 94 — 5 NE Rs +$ Ra — 3 i+ u | oosa 
7 1 7 E 1 à 3 

1 ra, cts Ev el,—2 3 60580 taht 

are S p 

zi en lo + q Kaa — 76 Li cos 9 
Vu = 2 ET cos] 
y. 2 4 1 
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Eier cilio nsa seo prege 
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da= E loq Kia — 2 Pit qi Kin 20080 
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In the formulae above are used the standard beta decay functions L;, M;, NL 
Q;, R; (Siegbahn 1955), the functions jg , M ; N 4 and R; of Curtis and Lewis and 
other new combinations of the electron radial wave functions. These new combinations 
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are defined in the Appendix I and some of them have been defined in the author's 
Z : 
previous paper. Limiting values for z4 «€ 1, po < l, but = ~ 1 (where p is the 
p 


electron momentum and o — the nuclear radius) have also been given. It is thus an 
easy task to write down the asymptotic form of our formulae’). 

In order to obtain a better review of the formidably long formulae it may be 
appropriate to start from the selection rules for different matrix elements. The results 
can be summarized in the following table: 


Type of transition Non-vanishing contribution from 
AJ = 0 (not 0 — 0) š all coefficients (except of Apg for 4 — 4 transitions) 
0—0 Án form =o 
AJ = 1 (not 0 —— 1) all Vp; and Ipis 4,5 and Apg 
0<> 1 all Vz; and Iki Aps 
AJ — 2 | Ars 


For a specific transition this short table enables us to consider only those terms 
in our formulae which do not vanish. 


The author is much indebted to Prof J. Werle for his interest in this work and 
for helpful comments. 


Appendix I. 


New B-decay functions. 


The following new combinations of Coulomb radial wave functions are introduced: 
! = [2p2 -1. 9-2. Sone 
Q6 = [2p? Fy (Z, E -0 as sin (4,— 4.) > aZ |> + 3 š +Ë 
aZ 


= Br* F (Z, Dt (fa — ae) sin (Ai A.) > ale ^ 
= Be RU BIS er fat ag) osli 4) > 3 Å 


2) A short discussion of the applicability of our approximation method and the comparison with 
the method used by Berestetsky et al. (1958a) is given in Appendix III. 
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Kia = [2p* Fy (Z, EN ose [f fy cos (4 ,— Ay) + 


1 
+ 818-2 cos (A, — 4.,)] > er 


K 12 = [2p? Fo (Z, E) Y OG [fi fr cos (d — 4») Sr 
2 
— 818-2 cos (44 — A_,)] > — = (2 is i 
Kj = Dp* Fy (Z E] - * [f., fy sin (A, — 4) + 
š aZ 
en (A, 4. > — AE PS 
K;-1,2 = [2p? Fy (Z, E)| 1 - 07? [f.1 fy sin (A, — A) — 
: lp 
— 18-2 sin (4, — A_,)] > — il 3E 
D Dorn ZEN o-* lift; gi. cos (4, — 4-4) + 
1 på 


81 fa cos (A, — Aj] M EYE ar 


aci [2p? F, (Z, EN e [fa g-a cos (4.,— 4.4) — 
1 1 på 
— & fa cos (4, — 49] >— | 3 pš + gp 
et SET [EY as o [fy Bro sim (42 — 4.) + 
+ 81 fa sin (A, — A>)] SE E+ n 
Ri i-a = [2p? Fy (Z, E)| 1: 0? [f.1 8-2 sin (4.5 — 4.9 — 


afin (4 — Ag) > — aZ (4 Å ++ 


The arrow indicates the limiting values of these functions for 2 «€ 1 and 
p 


po < 1 (but £ ~ 1). Here o denotes the nuclear radius and £ = LE The functions 


f, and g, are the radial wave functions of Rose normalized to. one particle per unit 
sphere. The Coulomb phase shifts are defined as in the paper of Rose et al. (1952). 
It may be noticed, for future application, that the definitions of the K functions can 


be written in concise form: 
= (2p? Fy? 0 DID ff sin (4; — 4) + 8-48, sin (4.,— 4) 
| = (2p3 Fg tof O [ y gj sin (4; — 4) + g_;f_;sin (4.; — 4.2] 


Ki 


Ki; 
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One always has to take | i| < |j |. By changing the sine functions in the defini- 
and K#,.1(j) 


. . ` .` 3 ` . + 
tions above into cosine functions we obtain the functions Kz; 


is given by 
: LEE 
LG) = UI 3 [sign G) — 1] 


The K functions together with the other f decay functions mentioned in this paper 
(the P’, Q', N”, R” functions being generalized for arbitrary indices) exhaust the 
whole variety of bilinear forms which can be formed by means of f; e'i and g; ed, 


Appendix II 


The bilinear forms of the F, functions. 


In this paper use was made of Deloff's projection operator A(7,p) which trans- 
forms a plane wave spinor into an exact solution of the Dirac equation in a Coulomb 
field. This operator was given in Deloff's paper (1958) with the accuracy to terms 
linear in (pr). This accuracy is sufficient in our calculations of first forbidden beta 
decay. We take this operator in the following form: 


AG D = ME, + BF; + à (GF) Fa + BF E GET) 
j=0 
where the first terms of the expansion are used: 


= a IE ge 43 » Á 24 - 
2p > VE—1 


PM mue b here 
2p*r | VE +1 | 


yE--1 ren 


IT Spr 

This operator and the bilinear forms of the Fj, functions arising from squaring the 
matrix elements were used in the calculations in some preceding papers of the author. 

It is a general feature that in deriving the formulae for n-times forbidden transi- - 
tions all F}; functions with the indices i <n must be taken into account while for 
n-times forbidden unique transitions the two functions Fy, and Fj, are not necessary. 
The first sixteen bilinear forms to be presented below have thus appeared already in 
the calculations of first forbidden unique transitions (Kuchowicz 1959) while the 
reemaining twenty expressions are characteristic for first forbidden non-unique transi- 
tions. All these bilinear forms are expressed in terms of the B-decay functions L; etc. 


[E Fy Z E)| (| Fiol? + | Fool? + E+ 2Re (Fio Foo)} = P, 
[E F (Z, E) 1 {E (| Fio iH | Fal?) + 2Re (Fio F,o} = L 
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[E R, (Z EI Fal? —[Faol?} = = Ly 
[E Fo (Z, E)|-1: 21m (Eis F; o) = = B 
[E Fo (Z, E) 1 {| FP + | Fs? + E-2Re (Fii F3) = 3 Pi 


[E Fy (Z, E) {E (| Fiz 


* 9 
2e ELS SE 2Re y Fay} = = pe 


0 
[E Fy (Z, E) 1 {| Fal? — [Fo 1l?} E = p 


* JT 

[E Fy (Z, E)! - 21m (Fi F1) = p Pi 
* * 6 5 

[E Fy (Z, ET - 2Re (Fii Fo = Foi Fso) = pp KE 
* 3 x, 6 5 

[E F, (Z, E) - 21m (Fii Fio = Fos F9) == på Kå,-2 


* * 6 = 
[E Fo (Z, E)]! - ZRe (F51 F0 — Fia Foo) = — RO 
* * 6 T. 
[E Fy (Z, E)|-* - 21m (Fo Fio = Fia Fao) = pi Ka 
* * * * 6 pom 
[E Fy (Z, E) : 2Re [Fi Fo "y Foi F2 E (Foi Fio HN Fi F»9)] = P Ka, 
* 3 * * * 6 — 
[E F, (Z, E) - ZIm [Fi F 1,0 + Foi Fz EE (Fai Fio XEM F,o)] == = Saa 
* * 6 
[E F, (Z, FIG « 2Re [E (Fi Fo Js Fri F0) "n Foi Fio im Fi Faol ~ 1 p Ki 
* * * * 6 
[E Fy (Z, E)! : 21m [E (Fii Fio + Foy Fs) + F31 Fio + Fin Fo] = + p Kia 
[E F, (Z, E) + {| Foil? + | Faal? + E: 2Re (F5, Fad} = — Qo — Pi + 2K 12 


1 , , gz. 
x = p (Mo +L, + 2K 1,2) 


[E Fo (Z, E) ^: {| Fal? SE | Fia 


[E Fy (Z, E]? - CE (Fal) + T 2) + 2Re (F4 F43) = Mo + L, + 2K;_ 10 


* 1 , , = 
[E F, (Z, E)]' : 2Im (F4, Fz) = 5 (Qo + P, +2K , 3) 
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RN WS 2 ; J 
LE Fi (Z E)! -2Re (F31 Fio e 4,1 F, 0) = p (act Kj) 


* + * 3 2 " 
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* A* 2 n den 
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Appendix III 


On the approximation method for relativistic coulomb wave function 


If one writes down explicitly the asymptotic form of our polarization-correlation 
az 

formulae for — > 1 and po < 1 (which was omitted because of lack of space, but 
P 


can easily be performed by inserting the limiting values given in this paper and in the 
references cited), some discrepancies between. them and the formulae for 0— 0 
transitions of the paper of Berestetsky et al (19582) will be found. If the substitutions 
employed by Berestetsky et al (19583) are made in eq. (5.8) of that paper, the results 
for the case of a pure Coulomb field should agreee with ours. This, however, does not 
happen. Apart from the question of opposite signs for the coefficients 4,5 and Aya 
(which is caused by taking the y; — matrix with the opposite sign by definition) 
the form of our coefficients 4,3, 44,4, Ass differs from the corresponding expressions 
in the above mentioned paper. There is also a difference (in the term which contains p? 
and in the signs of some other terms) between the expression which can be formed 
from ®, of eq. (5.8) for the case of Å interaction and our expressions. While the sign 
difference can be attributed to a possible misprint, the difference of the 4,4 coefficients 
seems to result perhaps from a different computational procedure. It is thus necessary 
to make clear what is the difference between the approximate Coulomb wave functions, 
which have a different form in this paper and in the paper of Berestetsky. 

Berestetsky's electron wave function can be written in the form that enables 
comparison with the function in Deloff's notation: 


VW Au — [Fro + ip: 7 Fa + Boo tip T Fay) +i@ T (Fa BF as) 


where u is a plane wave spinor and the functions Fj; have the form: 


| x ER zd ; 
Foo 2 VE + 1 VE — 1 F, 1 V2E E B 

NAME IL. USER c e 
; 7 Fa, Pal E$ zu 
FT 5 À e AR 
Fos 2; JET EE 


The coefficients z; and fj; are defined in the work of Berestetsky et al. (1958a). For 
a pure Coulomb field they have been given in a preceding paper of Berestetsky et al. 
(1958b). Comparison of the expressions obtained with the corresponding expressions 
of Appendix II gives the following results. Deloff's operator Å which was used 
throughout this paper, and Berestetsky's operator A’ (constructed in a similar 
manner to the former) are not fully equivalent. While in the limit for aZ — 0 
Deloff's operator goes over into 


/ 1 +ip:r 
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for first forbidden transitions, the corresponding operator A’ goes over into a form 
which contains also the term (x:7) multiplied by an expression linear in pr. This 
second term seems to cause the mentioned difference between our asymptotic values 
and those given by Beresteisky et al. (1958). This difference is evident only in the 
interference terms between the relativistic and retarded matrix elements, where the 
products of pr are of interest. 

In Berestetsky's eq. (1.2) are retained all terms of first order in pr. For the plane 
wave spinor representation, however, there is no need to retain terms of the first order in 
the functions F4, and Fj, because these functions are multiplied by ix: T and act 
on the small components of the bispinor. In our formulation the whole product is 
thusof the order (pr), while the wave function of the formula (3) in the paper of 
Berestetsky et al. (1958b) implies that also terms of a higher order are retained. 
The question of the omission of higher order terms will be discussed in a forth- 
coming paper. 

Not only these terms account for the difference mentioned. It should be emphasi- 
zed that the form of the radial wave functions in the paper of Berestetsky at al. (1958a) 
is hardly to be compared with ours. Averages over the nuclear volume have been 
also made by Berestetsky. It may be noticed that by inserting the expressions for F; A 
into the formulae from Appendix II we obtain the formulae (5.8) of Berestetsky et al. 
This was tested for the example of the expression for Fç. 

Note added in proof. The formulae presented in this paper can be adapted 
to allowed beta decay (including relativistic terms) after following substitutions: 


¿> f š fox? XT 
fi> fy z BA, 
Cy Z C, 
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LA CONSTANTE DIELECTRIQUE DES LIQUIDES DIELECTRI- 
QUES DANS LES CHAMPS ELECTRIQUES INTENSES. APERCU 
THEORIQUE ET EXPERIMENTAL *) 


Par A. PIEKARA 


Institut de Physique de l'Académie Polonaise des Sciences; Université A. Mickiewicz, Poznan 
(Recu le 16 février 1959) 


On observe trois types de l'effet de saturation diélectrique dans les liquides polaires 
et dans leurs solutions dans les solvents non polaires. Le premier est celui de l'effet négatif 
(diminution de la constante diélectrique sous l'influence du champ électrique) conservant 
son caractére négatif pour toutes les concentrations des solutions. Ce type d'effet se ma- 
nifeste dans l'éther éthylique, le bromure d'éthyle, le o-nitroanisol etc. Le deuxiéme type 
de l'effet de saturation — c'est l'effet positif avec inversion (nitrobenzéne, o- et m-nitroto- 
luéne), qui devient négatif dans les solutions de concentrations suffisamment faibles. 
Finalement, le troisième type de saturation présente l'effet positif sans inversion, effet 
qui reste positif dans les solutions de toutes concentrations. Ce dernier effet apparaît 
dans certains liquides dont les molécules présentent la rotation interne des groupes polaires, 
comme par exemple dans le 1,2-dichloroéthane etc. Le premier effet, négatif pour toutes 
les concentrations, est dà à la saturation proprement dite, tandis que les deux effets 
positifs sont provoqués par les interactions intermoléculaires (p.ex. dans le nitrobenzéne) 
ou intramoléculaires (p. ex. dans le 1,2-dichloroéthane). Dans ces deux derniers cas, 
c'est l'énergie de couplage des pairs des dipóles voisins qui, dans un champ électrique 
extérieur, fait naître un effet concurrent à l'effet négatif de la saturation normale 
(„ciseaux dipolaires^ s’ouvrant et se refermant sous l'influence du champ électrique 
extérieur, donnant lieu statistiquement à l'effet de > 0). 


Un gas aux molécules polaires, placé dans un champ électrique intense, prése- 
nte l'effet de saturation diélectrique, effet qui se manifeste par une diminution 4e 
de la constante diélectrique. Celui-ci, calculé par Debye d'abord pour les gaz et ensuite 
pour les liquides (voir Debye 1925, 1935), fut observé pour la premiére fois par 
Herveg (1920) dans l'éther éthylique. Le tableau 1 donne un aperçu chronologique 
des recherches ultérieures. (Les traits indiquent que les mesures se sont avérées 
impossibles à cause de la conductibilité électrique trop élevée). 

En parcourant ce tableau, on remarquera que la diminution de la constante 
diélectrique du nitrobenzéne trouvée par Malsch (1928) ne s'est pas maintenue, et 

*) Communication presentée lors du Congrés Diélectrique de l'URSS à Moscou en Novembre 1958. 
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Tableau I. 
eile asss Ehr (x benzène | Eau, alcool 
Auteur | Méthode éthylique | Nitrobenzéne au, alcools 
Je Herweg s 1920 | des battements (audio) | Ae < 0 — | — 
J. Kautzsch ...--. 1928 | des battements (audio) Ae < 0 — | — 
Ja Malsch. . «asi 1928 | d’etincelles Aree Ae <0 
H. Gundermann .. 1930 | d’étincelles 71613 00 El Ae < 0 
A. Piekara, | 
B. Piekara ...... 1936 | de résonances Ae « 0 Ae > 0 = 
ee Ed a == 
A. Piekara, 
A. Lempicki .... 1939 | des battements (oscillogr.) Ae<0 Ae > 0 == 
A. Piekara, 
A. Chelkowski .. 1956 | des battements (oscillogr.) Ae < 0 Ae > 0 = 


ceci par le fait qu'en 1936 nous avons trouvé que, dans le nitrobenzène trés pur, l'effet 
Ae est positif, c'est-à-dire que ce liquide placé dans un champ électrique présente 
une augmentation de la constante diélectrique. En méme temps nous avons pu 
constater que, par suite de la conductibilité électrique, il y a augmentation de la 
température dans le nitrobenzène insuffisamment purifié, ce qui entraîne une diminu- 
tion de la constante diélectrique; toutefois, ce phénoméne n'a rien de commun avec 
l'effet de saturation diélectrique. Ensuite, l'effet de saturation positif (de > 0) 
a été confirmé et retrouvé par Å. Chelkowski (1958) dans quelques autres liquides. 

Fig. 1. montre la méthode expérimentale. Deux générateurs, dont l'un de mesure 
à fréquence variable et l'autre, stabilisé avec un quartz (f 1 MHz), et qui sert de 


CRYSTAL 
OSCILLATOR 


AUDIO 
PLIFIER 


Fig. 1. Diagramme du montage pour les mesures de Ae = £ — ë, (e; — constante diélectrique du liquide 
étudié remplissant le condensateur C,; € — constante diélectrique du méme liquide soumis à l'effet du 
champ électrique extérieur). 
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comparateur, donnent des battements de fréquence acoustique, qui, après détection 
et amplification, sont transmis à une paire de plaques d’un oscillographe cathodique. 
L'autre paire enregistre une tension de la même fréquence provenant d'un générateur 
RC. Sur l'écran de la lampe oscillographique on observe une ellipse stationnaire. Le 
condensateur C, contenant le liquide étudié fait partie du circuit d'un générateur 


MO 
— Oo 


o 5 10 75 20 cm 


Fig. 2. Condensateur à liquide. 


de mesure. Aprés la mise en circuit de la tension de la source extérieure, la fréquence . 
subit un changement, qui est compensé avec un condensateur de mesure spécial de 
petite capacité. Fig. 2. montre l'un des types des condensateurs destinés au liquide 
étudié. C'est un condensateur extrémement solide en acier inoxydable dont la plaque 
intérieure est enfermée dans un cylindre creux. Ce condensateur résiste aux déforma- 
tions et sa construction permet d'y introduire le liquide par distillation dans le vide. 

En ce qui concerne le mécanisme de l'effet de saturation diélectrique dans les 
liquides, c'est l'étude de la grandeur 4e dans les solutions des liquides polaires dans 
des solvants nonpolaires (pour lesquels Ae = 0) qui présente un intérêt tout spécial. 
Or, d'aprés nos connaissances actuelles, il y a trois effets de saturation diélectrique, 
notamment (Fig. 3): 1. un effet négatif, Je < 0, qui se maintient dans le liquide pur, 
ainsi que dans ses solutions de concentration arbitraire; un effet de ce type a été 
trouvé dans l'éther, le bromure d'éthyle, le o-nitroanisol et dans d'autres liquides; 
2. un effet positif, Ae > 0, dans le liquide pur, qui change de signe au-dessous d'une 
certaine concentration; nous le désignerons ,effet positif avec inversion“; 
c'est le cas du nitrobenzène, ainsi que du orto- et méta-nitrotoluéne; 3. un effet positif, 
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Ae > 0, qui reste positif dans les solutions du liquide dans un solvant non-polaire 

å å : M 
de concentration arbitraire; nous le désignerons „effet positif sans inversion ; 
un tel effet a été observé dans du 1,2-dichloroéthane ou 1,2-dibromoéthane. 
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Fig. 3. Les trois types de l'effet de saturation (A=) observés dans les solutions des liquides polaires dans 

un solvant non-polaire: 1. Effef négatif (éther, o-nitroanisol etc.), 2. Effet positif avec inversion (nitro- 

benzène, o- et m-nitrotoluène), 3. Effet positif sans inversion (1,2-dichloroéthane, 1,2-dibromoéthane etc.). 
Inteusité du champ électrique extérieur: E = 60 kV/cm 


Le mécanisme moléculaire de l'effet de saturation positif avec inversion a trouvé 
son explication en 1937 (Piekara 1937, 1939). Nous admettons que c'est la molécule 
la plus proche qui joue un róle spécialement important dans l'interaction moléculaire 
d'une molécule donnée avec le milieu. Dans un liquide tel que le nitrobenzène, dont 
la molécule est fortement polaire, des paires transitoires de molécules se forment. 
Leur énergie de couplage (W) présente un minimum dans la position antiparalléle 
des dipôles des deux molécules. Le champ électrique extérieur agit sur eux de manière 
à diminuer quelque peu, en moyenne, l'angle entre les dipóles, ce qui conduit à une 
augmentation du moment effectif de la paire, et, de ce fait, à un accroissement de la 
constante diélectrique du liquide. Si cet accroissement l'emporte sur la diminution: 
de la constante diélectrique due à l'orientation, on observe l'effet positif. 
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Cette hypothèse conduit à la formule suivante pour la valeur moyenne statistique 
de la composante du moment électrique de la molécule dans la direction du vecteur 
du champ: 


ua F u^ F3 
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F étant l'intensité du champ interne. Ici, les „facteurs de corrélation“ de la polarisa- 
tion (R,) et de la saturation (R,) s'expriment par les relations suivantes: 


een 
L 
R; = Spr rrt P (1 + 9L) (3) 


où L = L (y) est la fonction de Langevin de la variable 


W 


. l'énergie de couplage en unités sans dimensions. La manière dont dépend R, de y est 
représentée par la courbe de Fig. 4. On voit que pour la valeur de l'énergie de couplage 
y — 1,33 le facteur de corrélation de la saturation change de signe. Ceci justifie l'exi- 
stence de l'effet positif avec inversion, puisque la variable y croít avec la concentration. 


Fig. 4. Le facteur de corrélation R, en fonction de l'énergie de couplage des dipóles. 
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Une comparaison des résultats théoriques et expérimentaux montre qu'il y a un trés 
bon accord non-seulement en ce qui concerne l'effet de saturation, mais aussi pour 
la polarisation et des effets d'orientation moléculaire tels que la biréfringence électrique 
ou magnétique (Piekara 1947, 1950). 

S'il s’agit de l'effet de saturation diélectrique positive sans inversion, le signe de 
cet effet, ne dépendant pas de la concentration, est indépendant de la distance qui 
sépare les molécules en solution; ceci indique que nous avons ici un phénoméne 


CI 


| 
intramoléculaire. En effet, dans le 1,2-dichloroéthane Ic — CH, il y a deux groupes 
Cl 
polaires qui peuvent pivoter autour de la liaison C—C. Toutefois, ceci présente un 
cas de rotation génée, puisque à la configuration trans (dipóles antiparalléles) correspond 
un minimum de l'énergie potentielle. Donc, cette situation rassemble à celle qui 
existe dans les paires moléculaires dans le nitrobenzéne; toutefois, les dipóles couplés, 
ici, font partie d'une méme molécule, ce qui fait qu'aucun degré de dilution ne saurait 


influencer sérieusement l'énergie de couplage des moments C—CI ni conduire à un 
changement du signe de l'effet de saturation. 


Dans la série des dérivés l,n-dihalogénés des hydrocarbures saturés, la distance 
séparant les moments des liaisons C—Cl (ou C—Br) augmente avec le nombre d'ato- 
mes de carbone dans la chaine. Or, dans le but d'apprendre si une augmentation de 
la distance entre les liaisons C—Cl (ce qui équivaut å une diminution de l'énergie de 
couplage de leurs moments) induirait un changement de signe de l'effet de saturation 
pour donner l'effet négatif, A. Chelkowski (1959), travaillant dans notre Institut, 
mesura l'effet Ae dans la série des dérivés 1,n-dihalogénés de l'éthane, du propane, 
du hutane etc. La Fig. 5 montre les résultats obtenus par lui en fonction du nombre n 
d'atomes de carbone dans la molécule. On voit que, pour les valeurs consécutives 
de n, il résulte, pour Ae, une ligne en zigzag. Ce qui est spécialement caractéristique, 
c'est le changement du signe de l'effet Ae, qui devient négatif en passant de n = 2 
à n = 3, et redevient positif pour n = 4. Ce fait doit être attribué à la géométrie des 
liaisons C—C, qui forment un angle d'environ 110°. Ceci a pour conséquence que, dans 
une molécule à deux atomes de carbone, les moments des liaisons C—Cl peuvent 
pivoter l'un par rapport à l'autre en traçant deux cônes très obtus ayant le méme 
axe (Fig. 6). Leur configuration privilégiée, c'est la configuration antiparalléle; quand 
l'énergie de couplage est suffisamment grande, ceci provoque un changement du signe 
de R. (Fig. 4) et conduit à l'effet. positif de saturation. Dans une molécule à trois 
atomes de carbone la situation est tout autre (Fig. 7). Ici, les axes des deux cónes 
forment un angle de 110°; donc, celles-ci sont à peu prés perpendiculaires; les moments 
des liaisons C—CI n'y présentent qu'un couplage très faible, et la configuration dans 
laquelle ils sont presque antiparalléles se réalise avec une probabilité beaucoup moins 
importante que dans le cas précédent. Par suite, on observe l'effet de saturation négatif. 
Dans la molécule à quatre atomes de carbone, la situation se reproduit telle qu'elle 
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Fig. 5. Effet de saturation Ae dans les 1,-dichloroalkanes (®) et 1,n-dibromoalkanes (o) en fonction 
du nombre n d'atomes de carbone dans la chaîne; Intensité du champ électrique extérieur: 
E = 35 kV/cm; d'aprés A. Chelkowski (1959). 
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Fig. 6. La molécule du 1,2-dichloroéthane dans la configuration antiparalléle des moments électriques 
des liaisons C—CI. 


Fig. 7. La molécule du 1,3-dichloropropane dans deux configurations avec les moments. électriques 
i des liaisons C—Cl presque perpendiculaires. 
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était pour n — 2, toutefois avec participation des configurations possibles dans le cas 
n — 3. Par suite de la plus grande distance entre les liaisons C—CI (Fig. 8), le couplage 
antiparalléle est bien plus faible ici; donc, l'effet positif de saturation, lui aussi, est 
moins prononcé. En passant aux chaines plus longues, a mesure qu'augmente le nombre 
d’atomes de carbone, le couplage des moments des liaisons C— Cl s'affaiblit davantage 
et leur orientation mutuelle cesse d'influencer le signe de l'effet Ae, qui, par suite de la 
petite valeur de l'énergie de couplage, reste négatif. Donc, les deux courbes P et N 


Fig. 8. La molécule du 1,4-dichlorobutane dans deux configurations avec les moments C—CI anti- 
paralléles. 


correspondant aux nombres pairs et impairs d'atomes de carbone convergent au- 
dessous de l'axe des abscisses (Fig. 5). 

Le point de vue que nous venons de présenter sur „lallure en zigzag“ de l'effet 
de saturation en fonction du nombre n se trouve confirmé si l'on compare les résultats 
des mesures des moments dipolaires (4,14) des dérivés l,n-dihalogénés des hydro- 
earbures saturés avec les valeurs théoriques de ces moments (uJ) calculées à partir 
de l'hypothése de rotation libre des liaisons C—Br (C. P. Smyth, 1955). La valeur 
de la difference Meale — Hobs peut servir de mesure du couplage entre les deux liaisons: 


1,2-dibromoéthane 


1,3-dibromopropane 0,3 
1,4-dibromobutane 0,7 
1,5-dibromopentane 0,4 


1,6-dibromohexane 0,3 
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On remarquera que, pour n — 3, le couplage entre les moments des liaisons C—Br 
est aussi petit que pour les grandes valeurs de n, mais que, pour les deux premières 
valeurs paires de n, ce couplage est exceptionnellement grand; ceci s’accorde de ma- 
nière qualitative avec l'allure de 4e en fonction de n (Fig. 5). - 


Pa 


CH, 


O 
Fig. 9. La molécule du o-nitroanisol; le radical // NCH, empèche le couplage entre les anneaux 
benzèniques de deux molécules voisines. 


Que c’est bien la rotation gênée autour de la liaison C—C qui joue le rôle décisif 
dans l'effet de saturation positif sans inversion (Fig. 3, courbe 3), c’est ce dont A. Chel- 
kowski a pu se convaincre en étudiant le comportement de Ae dans le 1,2-dichloro- 
éthyléne, dans lequel les deux groupes C—CI sont unis par une liaison double: HC— CH; 


el 
a ca 


or, il a constaté que Ae présente une valeur normale, donc négative. Ce point de vue 
est confirmé par les mesures cryoscopiques (Piekara, Chelkowski, Kielich, 1957). 
L'interaction des dipóles C—Br appartenant à des molécules distinctes ne joue pas 
de rôle décisif ici, puisque le monobromoéthane, C,HgBr, présente un effet négatif 
de Ae = —11.1075 pour E = 35 kV/cm. 

La rotation plus ou moins libre des radicaux dans la molécule ne favorise pas le 
couplage entre molécules dans les liquides. Pour s'en rendre compte, il suffit de con- 
sidérer le cas de l'éther, qui présente une polarisation diélectrique dont l'allure, en 
fonction de la concentration (p. ex. dans le benzéne), n'est point indicative de l'existence 
du couplage entre molécules; conformément à ceci, la saturation diélectrique 4e 
est négative. Plus éloquent encore est l'exemple fourni par une comparaison entre 
l'effet de saturation diélectrique dans le nitrobenzéne et celui présenté par le o-nitro- 
anisol. La molécule de ce dernier diffère de celle du nitrobenzène en ceci qu'elle com- 


porte un radical A CH, (Fig. 9) qui peut pivoter autour de la liaison C—O. Ce radical 
rend plus difficile le couplage des molécules aplaties pour former des paires transitoires, 
comme c'est le cas dans le nitrobenzène. Ceci expliquerait pourquoi l'effet de satura- 
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tion, qui est essentiellement positif dans le nitrobenzene, présente une valeur négative 


importante dans le o-nitroanisol (cf. Fig. 3). Le groupe polaire OT de par lui- 
méme, ne saurait produire d'effet positif, puisque son moment O—CH,, par ruite 
d'un effet stérique forme un angle différant peu de 90° avec le moment voisin 
C—NO,. 

Il existe actuellement plusieurs maniéres de traiter théoriquement la saturation 
diélectrique dans les liquides. Certaines de ces théories, toutefois, ne conduisent qu'à 
un effet de saturation négatif (Debye 1935, Frenkel 1945, Anselm 1944). La théorie 
la plus générale a été développée par A. D. Buckingham (1956) et par S. Kielich et 
moi (1958) en tenant compte de la déformation des molécules dans le champ électrique 
extérieur. Nous avons appliqué le méme formalisme aux molécules à rotation interne 
génée et nous avons montré qu'il rend possible une explication quantitative du 
phénoméne (Piekara et Kielich 1959). 
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A conception of a six-dimensional world enables a geometrical interpretation of 
the electric charge, charge conjugation, gauge transformations, and of the electromagnetic 
field. The same conception explains satisfactorily the isospin, its connection with the 
electric charge, and the ps-scalar character of nuclear forces. Several qualitative and some 
quantitative properties of strongly interacting particles (strangeness, rest masses, etc.) 
find an intuitive explanation within this geometrical framework. 


1. Introduction 


In the contemporary physical theories there exist two heterogeneous types of 
conservation laws: the conservation of energy, momentum, and angular momentum 
possess a geometrical background (homogeneity and isotropy of space-time) whereas 
the conservation of charge or of the baryon number possess none. This may be regarded 
as a methodological defect of modern theory and it may be used as a starting point 
for a more satisfactory development. 

The electric charge possesses a similarity to an angular momentum: the operator 
of an angular momentum component is proportional to the operator of infinitesimal 
rotation in a plane whereas the operator of the electric charge is proportional to the 
operator of an infinitesimal phase transformation. It is suggestive to regard this analogy 
as not being purely formal and to ascribe to the phase a meaning of an angle in a two- 
dimensional space, to be called „electric space“. In this way the electric charge ac- 
quires the meaning of a spin, whereby gauge transformations and charge conjugation 
are interpreted as rotations and inversion in the electric space. 

After ascribing physical reality to the electric space, the question arises as to 
the mutual relation between this space and the Minkowski space. The most straightfor- 
ward assumption will be to regard them as two subspaces of a six-dimensional metric 
space. We shall see that this assumption leads to a geometrical interpretation of the 
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electromagnetic field. The electromagnetic field will be amalgamated with geometry 
in a very similar sense as it has been done with the gravitational field: this produces 
also a curvature but in a six-dimensional manifold. The same conception of a six- 
dimensional metric space will enable us to explain the isospin and the connection 
of its *third* component with the electric charge. Several qualitative and some quanti- 
tative properties of strongly interacting particles (ps-scalar character of nuclear for- 
ces, strangeness, rest masses) find an explanation within the six-dimensional theory. 
The advantages of the hypothesis of a six-dimensional space are so conspicuous that 
they counterbalance the conceptual difficulties connected with the very fact that the 
fifth and sixth dimensions are not directly perceptible. 

The idea of interpreting geometrically the electromagnetic field within a six- 
dimensional framework is not new. Several authors worked on these linest). Although 
the very idea of accomodating the electromagnetic field in a six-dimensional manifold 
seems correct, some details of the former attempts do not seem satisfactory. In parti- 
cular, Podolañski assumed a ps-euclidean electric plane and investigated a two-para- 
meter gauge group interpretable as the group of translations in this plane. This leads 
to two "electromagnetic fields, one of them producing unphysical effects (owing 
to the indefinite character of the electric plane). Our presentation differs essentially 
from that of Podolanski. We investigate a six-dimensional space with a signature 
+ + + — — — so that the two electric dimensions are time-like, and limit ourselves 
to a one-parameter group of gauge transformations interpretable as rotations around 
a distinguished point. In this way we avoid Podolaüski's non-physical “second electro- 
magnetic field“. 


2. The distinguished role of the fifth and sixth dimensions 


By assuming a many-dimensional space, we are faced with the difficulty of explain- 
ing the obviously different role played by the additional dimensions as compared 
with the ordinary four dimensions. We can offer two alternative ,explanations*. The 
first of them is as follows: let us point out that also the time plays an obviously distin- 
guished role in comparison with the three space-like dimensions: we can perform arbi- 
trary translations in the three space-like dimensions but cannot do the same in the 
direction of the time except in our imagination. In reality we are “prisoners of our 
actual present“ and the actual “now“ is perhaps not a point but an infinitesimal ele- 
ment dx‘ on the time axis. Let us assume a “generalized now“ consisting in a three- 
dimensional element dx#dx5dx6 in the time-like subspace of the six-dimensional world. 
This element changes constantly and unavoidably in a certain direction in the time- 
like subspace, so that the world of our experience (present, past, and future) is not 
a four-dimensional manifold of points but a four dimensional manifold of two-dimen- 
sional infinitesimal elements (a world-tube). The Minkowski subspace can be characte- 


1) See Podolañski (1950), where also references to earlier works on this solied may be found. 
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rized by x9 = x$ = 0, but the world of our experience reaches beyond the Minkowski 
space and includes its infinitesimal neighbourhood in the directions of x? and x$. 
The acceptable transformations are those mapping the immediate neighbourhood 
of x9 — x9 — 0 upon itself. Thus, we can perform infinitesimal translations in the 
xŠ and x directions but macroscopic translations in these directions are not admitted 
since they would lead us beyond the world of our experience. Thus we admit as a fun- 
damental postulate of physics: physical laws have to be covariant only under trans- 
formations mapping the immediate neighbourhood of x9? — x9 — 0 upon itself. 

Our inability of getting beyond the infinitesimal neighbourhood of x5 = xë = 0 
is to be regarded as similarly inexplicable but similarly true and elementary as the 
fact that we cannot change arbitrarily our position in time, e.g. cannot come back to 
any past moment of our life. 

The alternative explanation runs as follows: We can make the geometry itself 
responsible for the distinguished role played by the Minkowski subspace. This can 
be achieved most simply — without going beyond the traditional framework of Rie- 
mannian geometry — by assuming a special form of the metric tensor that distin- 
guishes suitably four among the six dimensions. For example, we can assume that the 
two-dimensional electric subspace is not a plane but a surface of revolution. For 
definitness let us assume an ellipsoid of revolution. This surface possesses two 
distinguihed points P and P' lying on the symmetry axis. The six-dimensional space is 
to be understood as a four-dimensional manifold of such ellipsoids so that the four- 
-dimensional manifold of points P is to be interpreted as “the Minkowski subspace“?. 

The ellipsoid world-model admits the group of general transformations in six 
dimensions but, owing to the anisotropic properties of this model, the laws of physics 
have to be covariant only under the following subgroups: all transformations in the 
Minkowski subspace and rotations about the symmetry axes of the ellipsoids. 

Of course, such assumptions concerning the structure of the world necessitate 
in turn a deeper justification, but at any stage of every new theory we have to introduce 
a priori assumptions that cannot be justified otherwise as by their expediency in 
explaining and ordering the experimental facts. 


3. Electric charge and. charge conjugation 2 


The line element for the six-dimensional ellipsoid world-model with a signa- 


ture +++ — — — is 
ds? = (da)? + (da)? + (dx)? — (dx)? + gss(da5)2 + 2gggdx9da* + ga (dx5)2 (1) 


2) Such a model admits the existence of another „antipodic“ Minkowski space. This is an interesting 
possibility; the roles of particles and antiparticles could be reversed in the antipodic world consisting 
_ of points P’. But we can also avoid the existence of a second Minkowski world by assuming a model 
consisting of one sheet of a two-sheet hyperboloid. The hyberboloid can also degenerate into a half-cone. 
The vertices of the half-cones would be points in the Minkowski world. Our results will prove to be quite 
independent of the special model of the electric surface. 
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where o? = (x9)? + (x9)? and a, b are the parameters of the ellipsoid. The Minkowski 
subspace is defined by 


x5 = xŠ = 0. (3) 


In view of the structural properties of such a space the laws of physics do not 
need to be covariant under arbitrary transformations in all of the six dimensions but 
only under arbitrary transformations in the Minkowski subspace and under rotations 
around the symmetry axis of the ellipsoid. This last group is identical with the group 
of gauge transformations and accounts for the conservation of the electric charge. 
But we can also perform different rotations around the symmetry axes of the ellipsoids 
attached to different points in the Minkowski subspace. In other words the angle of 
rotation f can be made dependent upon the point in the Minkowski subspace. This 
group of transformations is nothing else but the group of gauge transformations of 
the second kind. Physical laws have to be covariant also under this group of non-linear 
transformations 


x5 — x5 cos f(x) + x8 sin f (x*), x5 — — x’ sin f (xF) + xŠ cosf(x*) (A) 


A charged field is to be interpreted as a six-vector whose only non-vanishing compo- 
nents are 9° and gê. The usual complex notation means 


- = (5 + ig) = die (på — ip’), (5) 


because the transformation (4) induces 
p> get, gt > te. (6) 
It is seen from (5) that the charge conjugation 


BERNIE (7) 
is to be interpreted as an inversion of the axis 6. Hereby is the sixth axis by no means 


privileged as compared with the fifth one, since we can perform also a gauge transfor- 
mation consisting in a rotation through the angle z. 


The charge conjugation assumes the simple form (7) not only in the case of charged 
bosons but also in the case of charged fermions provided a Majorana representation 
of the y, matrices is used. Therefore the four-component complex spinors can be 
also interpreted as vectors with respect to reflections and rotations around the sym- 
metry axis in the electric subspace. In both cases (of bosons or fermions) the electric 
charge is a unit spin within the electric surface). 
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4. À geometrization of the electromagnetic field 


The theory of gravitation was prepared by the following steps: (i) First of all 
there was the Lorentz group. The equivalence of all coordinate systems going over 
one into the other by means of translations and rotations accounts for the conservation 
of energy-momentum and angular momentum. (ii) There was also a freedom of per- 
forming general non-linear transformations whereby a variable metric tensor was 
introduced. However, this tensor was subjected to restrictions guaranteing the flatness 
of the space. In this case g, described only fictitious gravitational potentials. (iii) The 
introduction of true gravitation consisted in dropping those restrictions upon gy 
whereby flat space turned over into a riemannian space. 

The above mentioned three steps can be closely imitated in the case of electro- 
magnetics: (1) First of all we have the group of gauge transformations of the first kind 
(rotations in the subspace (5,6) around the symmetry axis of the ellipsoid). This accounts 
for the conservation of charge. (ii) There is also a freedom of performing gauge trans- 
formations of the second kind. From the point of view of the six-dimensional space 
this means a non-linear transformation whereby?) the g,,, become functions of x* through 
the intermediary of the variable phase f and its gradient 


LE: (8) 
It will be convenient to introduce polar coordinates 
HP 0 COR CT rp sini (9) 


whereby (1) assumes the form 


ds? = (dx)? + (dx?) + (da?)* — (dx)? + g, (do)? + Eaa (dx? (1) 
with 
mr). Baa = — Q’. 2) 
In these coordinates the gauge transformation of the second kind is simply 
«>a+f(e) (10) 
so that (1) becomes 
ds? = dy, gx* dx! + gee (de)? + Exa (da)? + 28ra dx* da (11) 
. with 
£y = 9g — O Safi Era = —0* fk (12) 


3) Latin indices run from 1 to 4, Greek indices from 1 to 6. 


E 
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where 
dy = 1 for k= DEA DS da = = l and 6,,= 0 for ki. (12) 


Boo and Zaa are the same as before. 

(iii) Closely imitating the procedure of general relativity we drop the restriction (8), 
i.e. assume an arbitrary four vector field f, (x, ... x$). This means an encroachment 
upon the geometry of the six-dimensional world. A straightforward but lengthy compu- 
tation shows that the three-indices Christoffel symbols are 


Tis = fr» Tis = —fr> Tes = Ipe = Ta = 0, (13) 


at the points o = 0 defining the Minkowski subspace. Therefore the covariant deriva- 
tives of a six-vector whose only non-vanishing components are g> and 9$ become 


V = 8, 9? + LP Pi = 8, 6 — f, p. (14) 
Going over to the usual complex notation (5) for the electrically charged fields we get 
ga = (8,— ifj) P» V = (8, + ifa) 9. (15) 


This provides a geometrical interpretation of the well-known expressions 9, if, 
occurring in the theories of charged fields. These expressions mean nothing else but 
covariant derivatives of electric *semivectors* (5) with respect to the coordinates a* 
of the Minkowski subspace taken at the points x? — x9 — 0 belonging to the Minkowski 
subspace. À 

Let us check the commutability of such covariant derivatives. It can be easily 
shown that 


Par — Par = — i (91 fy — 2, fp) = — i Fy. (16) 


In this way it is demonstrated that an electromagnetic field Fj, introduces a curva- 
ture into our six-dimensional manifold (quite apart from the “primary“ curvature 
of the ellipsoid). This yields a geometrical interpretation of the electromagnetic field. 
By inspection of the formulae (15) and (16) it is clear that our results are inde- 
pendent of the particular ellipsoid model of the electric subspace. Any surface with 
rotational symmetry (to account for the gauge invariance) and with a distinguished 
point (to account for the privileged role of the points forming the Minkowski sub- 
space) will yield the same interpretation of the electromagnetic field F,, and its po- 
tentials f,.4) 
EE ne à Ia EE Der ME erbe ee ee 


*) Note added in proof In this paper the electromagnetic field is treated as a given, external 
field. In a sübseqüent paper (tobe püblished in Nuclear Physics) Maxwell equations are derived 
from a lagragean that is identical with the curvature. The ellipsoid model is to be replaced by 
a new model where gy, and g,, ara singular like Q ^. The parameter a is closely connected with 
the Sommerfeld constant. 
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9. An interpretation of the isospin and the ps-scalar character of pions 


In this section it will be shown that the hypothesis of a six-dimensional world is 
also useful for an explanation of the most important features of nuclear interactions: 
the meaning of isospin and the ps-scalar character of pions. 

Let us start with the traditional equation for the nucleon coupled to the pion field, 


4 + — 
(2 Ve 94 + ig yg v D) y = 0. (17) 


The mass term has been dropped for the sake of simplicity as the nucleonic mass is 
probably due to an interaction, so that the mass term does not need to appear in a fun- 
damental equation. The nucleonic field quantity is endowed with two indices y,,, 
where x — 1, ... 4 is a spinor index with respect to the Lorentz group and p = 1, 2 is 
the isospin index. 

Let us rewrite (17) in a different form taking advantage of the following eight- 
dimensional notation 


Vu = P, Pa = Pa Vai = Py Vai = V, (18) 
Piz = D Vas = Po Var = V^, Vas = Pa 


In this notation (17) assumes the interesting form 


55, iX, PE 0 š (19) 


u=5 
where 


p=01, p=, g = da. (20) 


The matrices I", (Greek indices run from 1 to 7, latin from 1 to 4) are Dirac matrices 
with eight rows and eight columns 


n- [^ fn | “ne Í, cu 2) 
Yk Y5 Vys 


m =P See} =— in Tg. Ie c (22) 


They satisfy the usual relations 
IEEE = 205 (23) 


"The new form (19) of the Dirac equation satisfied by the nucleon reveals a deeper 
sense of the isospin and its connection with the electric charge. It throws also a new 
light upon the question why pions are ps-scalar. 

Let us recall that; eight-component spinors are geometrical objects fitting most 
naturally into a six- or a seven-dimensional space (Cartan 1938). Thus, the appearance 
of the seven matrices I7, in (19) and the simplicity and elegance of the interaction term 
appearing in (19) in Oi EROS to the clumsy form (17) suggest that the physical 
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space is not four but six or seven-dimensional. We shall limit ourselves to a six-dimen- 
sional viewpoint. 

Assume that the Minkowski space (l,...4) is embedded in a six-dimensional 
space (1, ... 6). The matrices I's and I, are connected with the fifth and sixth dimen- 
sions in a similar way as T}, ... I", with the usual four dimensions, whereas I’, defined 
by (22) plays the same role in the six-dimensional world as y; does in the four-dimen- 
sional world. 

The pion field is a vector and the nucleonic field is a spinor in a three-dimensio- 
nal “space“ formed by two genuine dimensions 5 and 6 and by an “imitation of a di- 
mension* 7 brought about by the algebraic fact that if we have six matrices satisfying 
(23) then there must exist also a seventh matrix of the same type (the product of the 
remaining six). Since not only I, but also J’; and I, can be regarded as products 
of the remaining six matrices, neither of them is privileged which accounts for the 
appearance of an imitation of a three-dimensional space (5, 6, 7) in spite of the fact 
that the physical space is actually only six-dimensional. 

In view of this situation it is intelligible that some physical phenomena (e.g. nuclear 
interactions) can take advantage of the algebraic symmetry brought about by the 
existence of three subsidiary matrices J'y, Tẹ, 17 but some other phenomena (e.g. ele- 
ctromagnetic and weak interactions) can be sensitive to the fact that the space is actually 
only six-dimensional and that the seventh dimension is not genuine. In this way we 
can understand qualitatively the most enigmatic feature of “isospace“: why not all 
but only some laws of physics possess a symmetry with respect to ratations in “isospace“. 
Moreover, we see why the *third* component of isospin is connected with the electric 
charge: it means nothing else but a rotation in the electric subspace (5,6) in accordance 
with the formalism developed in the foregoing sections. 

We can understand also why pions must be ps-scalar. Since from a six-dimensional 
point of view the interaction involved in (19) constitutes the simplest and the most 
natural possibility in the case of a meson field forming a charged triplet, and a transi- 
tion back to the traditional notation (17) leads automatically to the appearance of 
a yg in the interaction term, we see that a ps-scalar interaction with a ps-scalar meson 
field means the simplest possibility, by far more natural than a scalar coupling with 
a scalar meson field. 


6. The problem of strong interactions involving heavy mesons and hyperons 


Heavy mesons and hyperons can be also fused into eight-component spinors. 
For example we can form an eight-component spinor out of products of the field 
quantities describing the hyperon Å and the meson K 

4K, = (AK) 42K, —(AK), AK = (AK), — A,K; = (AR), 
AK, = (AR), AK = (AK), A;K, = (AK), A,K, = (AK), (24) 


or a mixed spinor-vector out of products of the field quantities describing the hype- 
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ron >; and the meson K 
Bi (DR DER — OE SKE OK}, aK, = (DK), (25) 
Dik, = IN =D DUR, = (DR), aK, = (IR) 


where u = 5, 6, 7. Now we can form the following interaction energy densities 


8a 73598 54 «3 
INN, or ig 2 N Te (IR) (26) 


Of course, such expressions have to be hermitized and antisymmetrized. Replacing 
in (26) the nucleonic field quantity N, by the eight-component spinor &, describing 
the cascade hyperon we get the whole list of possible strong interactions involving 
K meson. 


7. The nucleonic number and the electric charge 
There exist two independent fourvectors that can be constructed with the aid 
of the nucleonic field quantity Y and satisfying the continuity equations 
DU TDR Push TR Io UO (27) 
where 


Ts = > (Ts Ts — Iy D. (28) 


Of course, both fourvectors (27) have to be properly antisymmetrized. The continuity 
equation satisfied by the first fourvector in (27) accounts for the conservation of the 
nucleonic number. À linear combination of both 


ju = 3 Pr, (l + Lise) V (29) 


is to be interpreted as the electric charge and current density, and the continuity 
equation satisfied by (29) accounts for the conservation of the, electric charge. 
We notice that the electric charge and current density involves the projection 
operator 
l 


y2 


which yields an interesting analogy with the case of weak interactions involving as 


(1 + ie) (29°) 


1 : 
well a projection operator V3 (1 + ys). Expressions involving (29^) are not covariant 


under inversions in the subspace (5.6). Therefore, the charge conjugation is no more 
given by the simple formula (7) valid in the Majorana representation, but is to be 
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generalized as follows 


P>YV*G, V*>GY, (30) 


where G means a matrix with eight rows and eight columns of the form 


1 , 
— (11) (30) 


and 1 means a unit matrix with four rows and four columns. 

In spite of different characters of transformation under inversions in (5,6), both 
nucleonic number and electric charge behave in the usual way under gauge transforma- 
tions and therefore can be interpreted as intrinsic angular momenta (unit spins) 
with respect to the group of rotations in the subspace (5.6). 


8. On the meaning of “strangeness“ 


The formula of Gell-Mann 


b s 
€ ° 31 
where q denotes the electric charge, b the baryon number, I, the “third“ component 
of isospin, and s denotes the strangeness, is purely phenomenological and lacks any 


deeper justification. It can be made intuitively clearer by dropping the one halves 
g=n+L+T (31) 


where n is the baryonic number for baryons, but also n — 1 for K-meson. Since 
the electric charge as well as the *third* component of isospin are interpretable as 
intrinsic angular momenta (spins) in the electric subspace, the same must be true for 
the quantity T, so that we can write the formula (31^) inthe form 


q = n + Is + Ts, (32) 


which means that the strongly interacting particles possess three types of intrinsic 
momenta in the electric subspace (5.6) and that the electric charge is a resultant of 
them. Only the resultant angular momentum in the electric subspace produces a cur- 
vature, which explains why there exists only one electromagnetic field instead of two 
or more. 

The quantities q and n on the one hand and the quantities I and T on the other 
seem to be closely related as both q and n possess only integral eigenvalues whereas 
Ig; nor Tę aussme (simultaneously) integral or half-integral values. Since spi- 
nors constitute a natural conception in spaces possessing at least three dimen- 
sions, we suspect that I;, and T; are projections of some spins existing in (at least) 
three dimensions, whereas q and n (assuming only integral eigenvalues) can be acco- 


A Six-Dimensional Riemannian Manifold 381 


modated in a two-dimensional space (5,6) without having recourse to further dimen- 
sions. We know already that isospin can be accomodated in a three-dimensional 
“space” (5. 6, 7), where the dimension 7 does not need to be regarded for genuine, 
but is “imitated“ by the algebraic fact of the existence of three supplementary matrices 
T's, 15, 17. We cannot accomodate the quantity T in the same “space“ since in this 
case we would have to do with a reducible representation and the particle could not 
be elementary. But within the six-dimensional space with a signature + + + — — — 
there exists a possibility of accomodating a new sort of spin: namely a spin in the 
time-like subspace (4, 5, 6). Since the groups of rotations in the subspaces (1, 2, 3) 
and (4, 5, 6) commute, the particle can be characterized simultaneously by two diffe- 
rent spins: a space-like spin and a time-like spin. The former is the ordinary spin 
whereas the latter, or rather its projection (that means an intrinsic momentum compo- 
nent in the subspace (5,6)) is to be identified with T;g. In this way we find an interpre- 
tation of the quantity Tg, as a component of a three-dimensional timelike spin. 

Regarding the quantity Tgg as a physical reality we see that the contribution 5/2 
to the nucleonic charge is rather accidental. As a matter of fact, there are three different 
contributions to the nucleonic charge: the isospin +1, the baryonic number 1 and 
the time-like spin which in the case of the nucleon contributes — 4. The new quantity 
T3, is to be regarded as a physically meaningful substitute for the phenomenological 
concept of strangeness. The rule of Gell-Mann can be now stated as follows: for strong 
and electromagnetic interactions AT;g = 0, for weak interactions AT;, = + 4, 
whereas for | AT;g | = 1 the interactions are very weak or none. 

In the next sections we shall find further arguments in favour of the usefulness 
of the concept of the time-like spin T5, replacing the old concept of strangeness. 


9. A systematization of strongly interacting particles 


Let us arrange mesons and baryons in two tables according to the following rules: 
(i) We arrange them from the left to the right according to the increasing values of Tg. 
In the case of equal values of T4, we put the particle further away from the centre 
of the table if the absolute value of the isospin of the particle is larger. (ii) Then we 
arrange the members of each charge multiplet again from left to right according to 
increasing values of I, In this way we get the following tables. 
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These tables reveal interesting regularities. To each particle there corresponds an 
antiparticle situated symmetrically with respect to the centre of the table. The centre 
of Table 1. is occupied by the neutral pion which is the only particle identical with 
its antiparticle. The most astonishing fact is that the particles arrange themselves in 
an increasing sequence of their mass values according to their distance from the centre! 
In Table 1. the centre is occupied by the lightest among all strongly interacting partic- 
les — the neutral pion. Next comes the charged pion which is already a little heavier. 
Then comes the charged K meson and finally the neutral K meson which again is 
a little heavier than its charged counterpart. The same regularity is to be seen in the 
table of baryons. Close to the centre is the proton which is the lightest particle among 
the baryons. Next comes the neutron, then the hyperon Å and afterwards the triplet Z 
again in a correct order of increasing masses: 2+ is the lightest and Z- the heaviest 
member of the triplet. The mass of £ is yet unknown, but the regularity exhibited 
by our tables encourages us to make a prediction: Æ will prove to be lighter than 
its charged counterpart. 

It is hard to believe that the above mentioned regularities are purely accidental. 
Most probably the rest mass is derivable from an eigenvalue problem. The main 
features of the mass spectra are determined by the quantum number | T5, |. The fine 
structure of the mass spectra is determined by the absolute value of the isospin 7 and 
the hyperfine structure (or rather a splitting of the mass terms similar to the Zeeman 
effect) by its projection Ig. Indeed, we see that the masses depend in the first instance 
upon | Tg | which is to be interpreted as the maximal projection of T. It is striking 
that only the lightest among the strongly interacting particles, the pion, is characterized 
by T = 0, so that its whole mass is to be interpreted as a fine structure effect brought 
about by the fact that its isospin is I = 1. The difference of masses of the X and A 
hyperons can be also interpreted as a fine structure effetc. The hyperfine splitting is 
characterized by the following rule: the mass increases a little if I is oriented in the 
same direction as T and vice versa. 


10. A further discussion of the time-like spin 


Adopting the interpretation of the quantity Tg, as a component of the time- 
like spin, i.e. of a spin in the time-like subspace (4, 5, 6), we notice a restriction upon 
its possible orientations in the case of hyperons. Since Tg, Tep Tas can be denoted 
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alternatively by T4, T, T, respectively, we see that the sense of the internal rotation 
defined by T;g = T, together with the direction of the time axis defines a left-handed 
helicity in the case of baryons and a right-handed helicity in the case of antibaryons. 
This property of baryons in the subspace (4, 5, 6) is quite similar to the helicity of the 


neutrino in ordinary space. Moreover, since in the case of Z and Z as well as in the 


case of Æ and Æ there are only two orientations of T, it is seen that the projection T, 
of the spin T upon the time axis assumes only the extreme values. Such restrictions 
upon the spin orientations are possible only in a ps-euclidean space and are connected 
with the singularity of the light cone. Indeed, the subspace (4, 5, 6) is a subspace 
of a six-dimensional ps-euclidean space. The light cone in our six-dimensional world 
is given by 


(dx1)2 + (da)? + (da)? — (da)? — (dad)? — (dx)? = 0 (33) 


in the immediate neighbourhood of the Minkowski subspace x? — x9 — 0. The restric- 
ted orientation of the time-like spin T along the time axis should be connected with 
the conception of a “motion“ with a “velocity of light“ in the subspace (4, 5, 6) in 
the direction of x* (hereby the roles of time-like and space-like dimensions appear 
interchanged). Such a “motion“ is characterized by (33) together with the restrictions 
dx = dx*.— 0 so that 


(dx)? + (dx?)? + (dx)? — (dx)? = 0. (33°) 


However, this last expression means nothing else but a motion with the velocity of 
light in the Minkowski space. Therefore, in order to justify the restricted orientations 
of the time-like spin T we have to look for an internal motion of baryons with the 
velocity of light. The existence of such a motion can be taken for granted since the 
velocity operator is given by the Dirac matrices multiplied by the velocity of light 
c& and assumes the eigenvalues +c. 

To conclude it can be stated that our conception of a six-dimensional world with 
a signature +++ — — — enablesus to interpret the enigmatic quantity T appearing 
in (31^) as the maximal projection of a spin defined in the time-like subspace (4, 5, 6). 
Its restricted orientation and its additivity finds an intuitive explanation with the help 
of the Schrödinger “Zitterbewegung“ executed by fermions according to the Dirac 
theory. 


11. The problem of rest masses of heavy particles 


It is usually believed that the rest mass is a consequence of interaction with a “self- 
field surrounding the particle. We share this opinion but we would like to point out 
that a final result of a satisfactory computation of self-masses can be very simple and 
intuitive in contradistinction to what seems to be suggested by the awkward methods 
of self-energy computations (perturbation method etc.) applied nowadays. 

In Section 9 we have seen that the rest masses arranged themselves quite auto- 
matically in an increasing sequence according to the distance of the particle from the 
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centre of our Tables 1 and 2. This shows that the mass depends upon the quantum 
numbers T = | Tyg |, I, and Ig, and therefore must have something to do with the 
groups of rotations in the time-like subspaces of the six-dimensional world. This seems 
philosophically satisfactory: Relativity has taught us that the total mass of the particle 


is not a constant but a function of its momentum m = | p? + må. Therefore it seems 
plausible that also the rest mass should be a function of some (internal) momenta. 
If so, we can regard the mass as a purely kinematical quantity. 

The increasing of mass together with T constitutes the main feature of the mass 
spectra. Let us neglect the fine structure of the mass terms dependent on the isospin 
and focus our attention upon the main effect. Our conjecture is that the mass is a simple 
function of T with an intuitive geometrical meaning. It is known from general relativity 
that the particle is not a point but is represented rather by a Schwarzshild sphere 
with a radius proportional to the mass. In view of the similarity of the two euclidean 
subspaces (1, 2, 3) and (4, 5, 6), it is plausible to assume that the particle cannot be 
a point in the time-like subspace either but is represented also by a sphere with a radius 
proportional to the mass 

nem. (34) 
where I is a constant with dimension of a length. In this way we get two „internal“ 
degrees of freedom: the angles g, on this sphere. In other words, the particle can 
be represented as a rotator situated in the time-like subspace of the six-dimensional 
world. | 

Let us denote the internal cartesian coordinates of the particle in the subspace 
(4, 5, 6) (or in the subspace formed by the proper time and the coordinates 5 and 6) by 


n=rcosgsind, ry=rsingsnd.  rg-—rcos 0. (35) 


A straightforward generalization of the Dirac equation for the case of a particle which 
is free in the usual space but is a rotator in the (r,,r,,r4) space is 


(vi p; + tt k) y = — ky, (36) 
where 
> . 9 JB TEMA dero TT 
res re pad (tre +), Su (37) 


The components of 7 are Pauli matrices in the space (r,, rs, rs). The term appearing 
to the right hand side of (36) is necessary to ensure the compatibility with the supple- 
mentary condition (34). As shown in detail in former publications (Rayski 1957a, b) 
equation (36) together with the supplementary condition (34) leads to the following 
eigenvalues for the masses 


> 


iU CU d at Joe 
"r= (ri) > T= >> (38) 


where T is an absolute value of a (half-integral) angular momentum in the space 
(rus Ta 73). According to Table 2 we assume that the case T = 1/2 applies to the nucleon 
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and the case T = 3/2 applies to the hyperon Æ, whence we find that I is the Compton 
wave length of the nucleon and that the mass of Æ is 


= y2 m, = 2598 m, (39) 


This value agrees well with the experimental value 2590 m,. 
For the case of integral T we start with the Klein-Gordon equation 


(p; p; + k k) v = y (40) 
which can be applied to particles with integral as well as halfintegral (ordinary) spin. 
Now, the mass eigenvalues are 


1 
mr = [ILE Does. iL .. (41) 
Assuming the same constant / as before we get 


4 
m, = y2 my, = 2185 m, (42) 


which agrees very well with the mass of the hyperon À which is characterized by T — 1 
according to Table 2. Of course, the above numerical results could be regarded as 
purely accidental, but the probability of such an accidental agreement is small. 

The above results constitute a confirmation of our conjecture that T' is to be 
interpreted as a spin a three-dimensional space involving the electric plane and oriented 
so that T;g is its extreme component. 

Å still more intuitive picture of the rest mass is obtained in a classical limit. 
Classical approximation is applicable to large values of the quantum numbers.Neglect- 
ing 1 or 1/2 in comparison with T both (38) and (41) go over into 


m = ua , (8) 


The fact that the mass is proportional to the square root of an angular momentum T 
can be understood in the following way: The particle structure is represented by a rota- 
tor in the timelike subspace. The simplest assumption concerning the structure of 
such a rotator is the reciprocal symmetry r = p in natural units c = Å = l = 1 where r 
is the radius and p is the momentum of a rotating point. Since r = m we get for the 
angular momentum T = rp = m? whence (43) follows. The particle is a three-dimen- 
sional rotator but oriented so that its rotation takes place only in the plane (5.6). 
'The sense of this rotation is such that together with the direction of the axis 4 it defines 
a left-handed (right-handed) screw for particles (antiparticles). 
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NUCLEAR SPIN-LATTICE RELAXATION TIME AND THE ELEC- 
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The intrinsic dielectric relaxation time 7 for polar dielectrics has not yet 
been measured. Dielectric loss measurements yield only information about the macros- 
copic dielectric relaxation time T, but they do not concern directly the intrinsic relaxa- 
tion time T. f 

Efforts have been made by different authors to calculate + as a function of T. 
These calculations are based on different internal field theories and they yield different 
results. In the Fig. 1 are shown the results of three theories (taken from: the paper 
of Powles 1953) in form of T/t versus ¢,/s,, (e, — static dielectric constant, Eo — 
dielectric constant in long wave infrared region). &, is taken to be 5.5, the value 
found for water (Collie et al. 1948). 

The aim of this paper is to obtain the intrinsic dielectric relaxation time by means 
of the Bloembergen-Purcell-Pound (1948) equation 7 = 37, (t, is the correlation 
time) for a given polar liquid at different temperatures (for the change of &,) and then 
to compare this values with the experimental values of T. 

In this way the ratio T/t for water at different temperatures was calculated and 
the results are shown as small-circles on the diagram (Fig. 1) The values of v, has 
been calculated by means of the Bloembergen-Purcell-Pound theory corrected by 
Kubo and Tomita (1954): 


1/ T, =; 1/ T, wn 1/ I 
1/T,, = 1.5 y4 h*v Jb = 7.03 x 1019 v, 


where T, is the nuclear magnetic spin- -lattice relaxation time, 1/ T, > 1/73, is the rota- 
tional and translational parts of 1/T,, y is the gyromagnetic ratio of the proton, 7, is 
correlation time, b in the proton-proton distance in the molecule. For water b — 
= 1.515 10-8 cm (Landolt-Börnstein 1951). The T, values are taken from previous 
measurements (Hennel et al. 1958). The interpolated values are shown in Table 1. 
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The estimated random error for interpolated 7, data is 2%. The values of T are taken 
from the measurements of Collie et al. (1948). 

A separate problem is how to devide the experimental values 1/T, into the rota- 
tional (due. to intermolecular spin-spin interaction) and the translational part (arising 
from interaction with protons situated in neighbouring molecules). According Blo- 
embergen et al. the ratio T,/Tj, is independent of the temperature. In this paper the 
value Tj/T,, = 0.478 is assumed so that the point for 20°C would lie on the Powles 
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Fig. 1. Relation between macroscopic T and intrinsic 7 dielectric relaxation times according to different 


internal field theories; a — Lorentz-Debye theory: T/r = (e; + 2)/(& + 2), b — Powles theory: - 


Tir = 3¢,/(2€, + £y) c — Onsager-Cole theory (Powles 1953). Experimental points o are obtained for 
water at different temperatures. 


curve (Fig. 1). But the change of this value (in sensible limits, of course) causes only 
a vertical shift of the experimental T/r curve without any greater influence on its slope. 
| The striking agreement of the experimental points and especially the slope with 
the Powles theoretical curve is worthy of consideration. It would be very interesting 
to extend the experimental data towards higher temperatures. 
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Table I. Interpolated values of T, for pure water 
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Acta Physica Polonica 18, 37 (1959). 
Equations of Motion of Rotating Bodies in General Relativity Theory by W. Tulezyjew. 


The terms (2.44) do not lead to correct equations of motion (2273 The 


correct Lagrangian contains second derivatives of 7” and €”. In order to avoid 
1 tek 2 


T a 1 
this complication we put S” =} SE 7° and 5' = å S" C^ in place of S" = 5" = 0. 
Now the terms connected with rotation in the post-Newtonian. piu of mo 
tion have the form: : 
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First ondio Beta Transitions for ihe Universal Fermi Pose 


by B. Kuchowicz 
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